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Abstract 

We present an intuitive, light front reduction scheme for the Bethe-Salpeter equation that consists of iterations 
of the (light front) energy dependent covariant equation followed by an instantaneous approximation. This scheme is 
equivalent to and motivates the form of a standard reduction that reproduces light front, time ordered perturbation 
theory for bound states. We use this three-dimensional light front formalism to compute electromagnetic form factors 
between bound states to show how higher, light front Fock components arise from the covariant Bethe-Salpeter wave 
function. Moreover, we show that non-wave function vertices, which arise after integration over the minus momenta, 
are completely removed in favor of higher Fock states when one correctly takes into account all contributions at a 
given order in perturbation theory. Such non-wave function vertices are shown only to occur in models in which 
the interactions are instantaneous. These models lack true higher Fock states. The normalization of the reduced 
wave function is derived from the covariant framework and related to nonvalence probabilities using familiar Fock 
space projection operators. Using a simple model, we obtain expressions for generalized parton distributions that 
are continuous. The nonvanishing of these distributions at the crossover between kinematic regimes (where the plus 
component of the struck quark's momentum is equal to the plus component of the momentum transfer) is tied to 
higher Fock components. Lastly we apply the light front reduction to timelike form factors and derive expressions for 
the generalized distribution amplitudes in this model. 

1 Introduction 

More than a half century ago, Dirac's paper on the forms of relativistic dynamics B introduced the front form Hamilto- 
nian approach. Applications to quantum mechanics and field theory were overlooked at the time due to the appearance 
of covariant perturbation theory, in which the special role played by time was supplanted by Lorentz invariance. Front 
form dynamics gradually resurfaced, first under the guise of the infinite momentum frame (in which covariant calcula- 
tions could be achieved more easily if the processes were viewed at the speed of light) and later as quantization on the 
light front plane (or null plane). 

The reemergence of front form dynamics was largely motivated by simplicity as well as physicality. The light front 
approach has the largest stability group of any Hamiltonian theory, that is light front dynamics maintains the smallest 
number of interaction dependent generators of the Poincare algebra. In highly relativistic calculations the utility is clear: 
Lorentz boosts on the light front are merely kinematical; whereas in the instant form, boosts depend on the dynamical 
evolution of the system. The physical motivation for hard scattering processes, such as deeply inelastic scattering, 
was the infinite momentum frame's utility in describing partonic subprocesses. Today this connection to light front 
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dynamics is transparent: hard scattering processes probe a light cone correlation of the fields. Not surprisingly, then, 
many perturbative QCD applications can be treated on the light front, see e.g. ||. Outside this realm, physics on the 
light cone has been extensively developed for non-perturbative QCD [Q as well as applied to nuclear physics || |(J . 

This paper concerns current matrix elements between bound states of two particles in the light front formalism. We 
approach the topic, however, from covariant perturbation theory. By projecting covariant quantities onto the light cone, 
one has a way of deriving light front amplitudes from field theory, circumventing the often subtle task of formulating 
the front form of dynamics ab initio. As demonstrated by the tremendous undertaking of Q , one can derive light front 
perturbation theory from covariant perturbation theory for scattering states, thereby demonstrating their equivalence — 
including the delicate issue of renormalization. As to the issue of light front bound states, a reduction scheme for 
the Bethe-Salpeter equation recently appeared || that produces a kernel calculated in light front perturbation theory. 
These authors demonstrate their reduction but they do not focus on its origin. We adopt this reduction scheme of 
the Bethe-Salpeter equation below, however, only after suitable intuitive motivation. For the purpose of simplicity, we 
consider only bound states of two scalar quarks interacting via the exchange of a massive scalar. The intuitive reduction 
scheme is quite similar to one found in however, we focus of the Bethe-Salpeter vertex function and its energy poles. 

Our main consideration is to extend the reduction to current matrix elements to investigate valence and nonvalence 
contributions in the light front reduced, Bethe-Salpeter formalism. Thus we calculate our model's form factor; moreover, 
we choose not to exploit a special reference frame in which we can neglect Z graphs. This enables us to completely 
investigate their contribution (as well as functional form) which has a variety of applications from timelike decays 
and timclike form factors to generalized parton distributions. These Z graph contributions have haunted light front 
dynamics since the nonvalence properties of the bound state are involved and thus valence wave function models cannot 
be utilized directly. A variety of solutions to this dilemma have been pursued. In JTo|, the nonvalence piece was treated 
as a new sort of wave function with an extra parameter determined by fitting to data. The approach allows for some 
phenomenological consistency, although the underlying field theoretic nature of the problem is ignored, e.g. Lorentz 
invariance is not maintained. The authors p"l| propose that the light front Bethe-Salpeter equation can be iterated to 
relate the nonvalence contribution to a valence wave function. In lieu of knowing the interaction which performs this 
crossing, they propose that the nonvalence vertex can be replaced by a constant. The reader is then riddled by their 
insistence that the constant vertex is consistent with data and therefore generally true, and later falsely assured that the 
vertex can uniquely be determined by Lorentz invariance (and hence depends only on one variable £, perhaps two i). 
The constant approximation and Lorentz invariance fitting are far from unique since the complete functional dependence 
(on £, t, x, k- 1 ) of the vertex is ignored (and hence true nonvalence correlations are wiped out). The procedure also results 
in conspicuously discontinuous distributions which one can trace directly to the constant vertex approximation. One 
could just as well choose a vertex with some arbitrary dependence on the neglected variables that forces continuity and 
includes an overall function to be determined from Lorentz invariance. This approach is a hybrid of Jl§ and that 
has not been pursued, until very recently Jl^] . 

We, too, tried a similar approach of iterating the reduced Bethe-Salpeter equation | fl3| in an attempt to rigorously 
generalize the results of [[lj} . We used the same simple Lagrangian model used below to express the crossed interactions 
and it seemed one could avoid explicitly using higher Fock space components. Our results, however, were not rigorously 
obtained in perturbation theory. We show below that contributions from non-wave function vertices are exactly canceled 
when one works completely to a given order in perturbation theory. Thus one cannot side step the issue of higher Fock 
components. 

We begin in section |^ by considering an intuitive reduction scheme for the Bethe-Salpeter equation. The method 
consists of a series of iterations of the covariant equation, followed by an instantaneous approximation for the vertex 
function, which allows for all integrals over minus momenta to be performed. The focus concerns the light front 
energy poles of the Bethe-Salpeter vertex, and further iteration before an instantaneous reduction allows for a better 
approximation to the poles. We use this scheme to derive the bound state equation to second order in the interaction, 
making the connection to light front, time ordered perturbation theory. We then use this intuitive scheme in section [|to 
motivate the formal reduction carried out by j8j. Once we have the formal scheme, we reverse the logic and demonstrate 
how the intuitive scheme arises. Discussion of the light front wave function's normalization is delayed until Appendix 
A. 

The formal reduction is then used to obtain current matrix elements between bound states in section ^|. We first 
formally construct the gauge invariant current and then apply it to form factors in perturbation theory. At leading 
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order we are confronted with non-wave function vertices, but in order to appeal to crossing symmetry we must go to 
higher order. Considering then the complete next to leading order expressions, we show the non-wave function vertex 
contributions are removed and replaced by higher Fock components. In Appendix B, we describe how crossing symmetry 
can be exploited in the case of an instantaneous interaction, such as QCDi + i |14) . 

The immediate application of our development for form factors (in a frame where the plus momentum transfer is 
non-zero) is to compute generalized parton distributions. We carry this out in section]^ where continuous behavior at the 
crossover is demonstrated. We also make explicit the relation to Fock component overlaps in this model. Connecting back 
to the intuitive scheme, we notice that the poles of the covariant vertex function generate the higher Fock components 
present in our expressions for form factors — and quite naturally the Bethe-Salpeter wave function retains all Fock 
components (in principle). This necessarily implies the same is certainly not true of the reduced light front wave 
function. Its lack of light front energy dependence restricts it to the valence sector alone. In Appendix C, we derive 
the higher Fock components for the model directly by using old fashioned, time ordered perturbation theory for the 
light front Hamiltonian. Another application is pursued in section [| where we derive expressions for the generalized 
distribution amplitude for this model, which is related to the timelike form factor via a sum rule. Such an application 
is interesting since no Fock space expansion in terms of two-body bound states is possible. Lastly we conclude briefly 
in section [?]. 



2 Intuitive reduction 

We start by writing down the covariant equation for the meson vertex function T. It satisfies a simple Bethe-Salpeter 
equation |15 , first schematically as operators (see Figure]]]): 



T = VGT, (1) 

where G is the full, renormalized, two-particle disconnected Green's function (which can be expressed as a product of 
two, renormalized, single-particle propagators). Here V is the interaction potential, which is the irreducible two-to- 
two scattering kernel. Let us denote the meson four-momentum by R, the quark's by fc, and the physical masses by 
R 2 = M 2 , fc 2 = to 2 . Since we take our quarks to be scalar particles, the renormalized, single particle propagator has a 
Klein-Gordon form: 

d ^ = (k 2 - to 2 )[1 + (fc 2 - m 2 )f(k 2 )} + ie' ^ 

where the residue is i at the physical mass pole and the function f(k 2 ) characterizes the renormalized, one-particle 
irreducible self interactions. We shall assume there are no poles (besides at the physical mass) in the propagator. This 
is consistent with model studies [|l6|. Since our analysis hinges on poles, we can neglect /(fc 2 ) in what follows by noting 
that to add it back in, we merely evaluate at the relevant poles (since all integrations in question will be calculated from 
residues). 

In the momentum representation, the Bethe-Salpeter equation then appears as 



n /,-./? i = / / AP V (k,p)*Bs(p,R), '••!) 



in which we have defined the Bethe-Salpeter wave function ^bs as 

* BS (k,R) = G(k,R)r(k,R), (4) 

with G{k,R) = d(k)d(R — fc). Since we will ultimately be concerned with soft physics probed in hard processes, we 
imagine the initial conditions of our system are specified on the hypersurface x + — 0. (We define the plus and minus 
components of any vector by x = -^(x ± x 3 ) .) Correlations between field operators evaluated at equal light cone time 
turn up in hard processes and thus this choice of initial surface is natural. In order to project out the initial conditions 
(wave functions, etc.) of our system, we must perform the integration over the Fourier conjugate to x + , namely fc - . 
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Figure 1: Diagrammatic representation of the Bethe-Salpeter equation. The blob represents the vertex function T. 



For instance, our concern is with the light front wave function defined as the projection of the covariant Bethe-Salpeter 
wave function onto x + = 0, 

1 p(k+,k ± ,R) = 2R+x(l - x) J -^^Bsik, R). (5) 

Looking at Eq. (Eh, in order to project the wave function exactly, we must know the analytic structure of the bound 
state vertex function. If the vertex function T(k, R) had no poles in k~, then our task would be simple: the light front 
projection of ^>bs would pick up contributions only from the poles of the propagator G(k,R). Next we observe from 
Eq. (||), that the k~ dependence of the interaction V(k,p) must give rise to the k~ poles of the vertex function T(k, R). 

If we imagine an interaction which is dominated by a piece instantaneous in light front time x + (i.e. independent of 
the light front energy k~), the vertex function will similarly be independent light front energy. Under these assumptions, 
the projection onto the initial surface can be carried out yielding a zeroth order approximation to the light front wave 
function. (The equal time version of this approximation was first carried out by Salpeter jl?]] but by specifically using 
the instantaneous Coulomb interaction.) Before worrying about how to correct this approximation, let us spell out the 
procedure and form of the zeroth order wave function. 

The vertex function T is replaced by its instantaneous approximation 7 which must be three-dimensional. Thus 
j(k, R) depends upon k 2 = 2 k+ k~ on — k ±2 = m 2 , R 2 — M 2 , and k ■ R, where 

k ■ R = k" on R+ + k+ R~ - k 1 - -R 1 - 

= ^(m 2 +xM 2 + k x2 Iel ), 

with x = k + /i? + and k re i = ^-xR 1 . We shall write this dependence as 7 — 7(2;, k^ei \M 2 ) and often for 
simplicity omit the relative subscript on the transverse momentum. The projection is then carried out simply as noted 
by integrating the propagator over k~, 

dk - m< R - k} - / 2k+ £ k _ k _ ^ (6) 



where 



^ a ^ on 2 k+ 

k b = R + (k — R) on — 2 (k+-_R+) ■ 



(7) 



The above integral is performed using the Cauchy integral theorem. For k + < both poles a and b lie in the upper-half 
complex k~ plane while for k + > R + both lie in the lower-half plane. In both cases, we can close the contour to avoid 
enclosing any poles and the integral vanishes. On the other hand, for < k + < R + closing the contour in the upper-half 
plane picks up 

J dk- d(k)d(R -k) = 2ni Res(k" 6 )6»[a;(l - x)] 

= ^ Dw(-T, k^^ \M 2 )9[X(1 - X)], (8) 
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where we have defined the Weinberg propagator as 



x(l — x) 

The zeroth order, light front wave function is then 

^(°) = i D w (x, k x rel |M 2 ) 7 (z, k x lel \M 2 )9[x(l - x)} (10) 

and hence is a function of only the relative momenta. 

To better approximate the wave function, we need to correct for using the instantaneous approximation. We can 
achieve this systematically by iterating the covariant Bethe-Salpeter equation for T. A better approximation to i\) would 
be iterating once and then replacing T with the instantaneous vertex 7. By retaining energy dependence in V , one 
is allowing for minus momentum poles in T, though we still must make an instantaneous approximation in order to 
carry out the integration. The next better approximation uses two iterations before using an instantaneous vertex 
replacement. Schematically the chain of successively better light front wave functions can be depicted as 

V^ (0) = J G 1 

?A (1) = / GVGj 



V-W = / (GV) 2 G 7 



Of course the exact light front wave function is ip 1 -- 00 ^ for which the poles of the vertex remain unaltered. At any 
level in the scheme, ipW will depend upon the j — 1 previous approximations. To derive a bound state integral equation 
for ipti', one merely replaces all ij)^ k \ k < j with rpv>. This is possible since in the equation for ipv>, a term with ij)^ k ' 
will be multiplied by j — k powers of V . Imagining that the interaction is weakly retarded, we can upgrade ip( k ' to ijj^' 
for free, since the difference introduces terms of order at least j + 1 in the interaction which are irrelevant: ipV' is valid 
only to order j.^j 

To illustrate this intuitive scheme, we desperately require an example. Indeed to say anything less than general, we 
must know the k~ dependence of the interaction. We therefore adopt a weakly coupled, one-boson exchange model for 
V (the so-called ladder approximation). Supposing the boson mass is [i and the coupling constant <?, we have 

V(k,p) - 



2(p+-k+)p--pV 
where the energy pole (with respect to p) of the interaction is 



(p — k) 2 — /i 2 + ie 

(11) 



-g 2 1 



_ , (p^k^+Zi 2 ie 



2(p+-k+) 2(p+-k+) 
With this model, we will work out two non-trivial orders of the reduction scheme 



1 In this way we are making a perturbative expansion for the reduced kernel and not counting powers of the interaction (which must be 
large GV ~ 1 to produce a bound state). This will be made clear for the ladder approximation to the kernel below. 
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Case 


y<0 


y<0 


0<y<l 


0<y<l 


y>i 


y>i 




y<x 


y>x 


y<x 


y>x 


y<x 


y>x 


Positive 


a b v 


a b 


b v 


b 


V 




Negative 




V 


a 


a v 


a b 


a b v 



Table 1: All possible combinations of signs of the poles' imaginary parts for the first order, light front wave function. 
The p poles are labeled as they appear in Eq. (nj 



2.1 First order interaction 

To first order in a weakly retarded interaction, we have 

^\xM) = 2R+x(l-x) J ^G(k,R)t J ^V(k,p)G(p,Rh(y,p^ rcl \M 2 ), (13) 

where we have defined y = p + / R + . Having inserted the instantaneous vertex 7 after one iteration of the Bethe-Salpeter 
equation, the integral over p~ can now be performed in addition to the k~ integral. To begin, first consider the p~ 
integration. The integrand has three poles (two from the propagator and one from the interaction) 

P a — P on — 2R+ y 

p- b = R- + ( P -R)- n - ^t^j (14) 

rT _ l r - 1 (p^-k-^+p 2 ie _j_ 

. P v — K 2(p+ - k+) 2BJ y-x ' 

which have imaginary parts that switch sign depending on the size of x and y. In Table |], we enumerate the possible 
combination of signs of the poles' imaginary parts. So for example, when y < and y < x, we close the contour in the 
lower-half plane and the integral vanishes, as is the case additionally when y > 1 and y > x. To simplify our work, we 
can assume that contributions to if) will require < x < lj^] Thus two regions are eliminated due to impossibilities: 
(y < 0) fl (y > x) = (j) and [y > 1) fl (y < x) = (f>. The integral is 



27ri%(l - y)} (% - x) Res(p- b ) - 0{x - y) Res(p-j) 



2R+x(l-x) j ^G(k,R) J _ y) ^( k ~> ^ P^H)^ (y, p' re i), (15) 



where we used 6[y(l - y)] D w (y, p 1 - \M 2 )~f(y, p x |M 2 ) = -iip(°\y, p ± \M 2 ) and defined 



F(k-, x, ; y, p^) = g' ( + ) ■ (16) 

2R+{y-x) yp b -p v p a -p V J 

Above we commented that tf>^ would depend upon the lower order approximation and this is indeed the case. To 
derive a bound state equation for the first order, light front wave function, we merely realize that changing ip(°' to ip^ 1 ' 
above results in corrections of second order in the retarded interaction. To derive this three-dimensional bound state 
equation, we need to perform the k~ integral. V retains k~ dependence through p - ^ and hence a renaming scheme is 
necessitated. Let us define 

_ (k^-p-^-tV I ie x 

\ K v a - 2R+{x-y) +P on+2fl+ (y-x)y n ? s 

i k _ _ o^-py+y R - , R) - u x-i V 7 ) 

l K «b — 2R+(x-y) ^ \l J rt /on T 2 _R+ ' 



2 



2 If we do not assume this, the additional residues with respect to p~ which we pick up in Table |lj turn out to evaluate to zero when 
integrated with respect to k . In this way, we need not assume the wave function tjiW is zero outside < X < 1; extra computation will 
demonstrate that it indeed vanishes. 
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Case 


x<0 

y<x 


x<0 
y>x 


0<x<l 
y<x 


0<x<l 

y>x 


X>1 

y<x 


X>1 

y>x 


Positive 


a b v a 


a b Vb 


b 


b v b 






Negative 






a v a 


a 


ab v a 


a b Vb 



Table 2: All possible combinations of signs of the poles' imaginary parts for the first order bound state equation. The 
k~ poles are labeled as they appear in Eqs. (0) and (|l7|). 




Figure 2: Diagrammatic representation of the one-boson exchange potential V appearing in Eq. (|20|). 



which renders 

^^> = ^(^ + ^)- 

Four poles now confront us in performing the second minus momentum integration: {k _ a ,k _ f,,k _ 1 , a ,k _ 1)6 }. All possible 
combinations of signs of their imaginary parts are listed in Table |^. Using Cauchy, we end up only with two residues to 
evaluate yielding the desired bound state equation 

V (1) (z,k\el) =D w (z,k\. cl |Af 2 ) / -^^^V(x,k\ e r,y,p\ el )^(y,p\ e} ), (19) 

J 2(2tt)^(1 - y) 

for the first order interaction V defined by 

V(x, k\ cl ; y, p\ ol ) = 9(x - y)V(k~ b , x, k x ; y, p x ) + % - x)V(k~ a ,x, k x ; y, p x ). (20) 

The equation takes the form of the Weinberg equation jl8| for the light front time ordered, one-boson exchange potential, 
since, for example 

^A.v ik -.,^-,y, P -r - R- - r- - fr^y' + / - (R-pK. (2D 

is the energy denominator corresponding to the intermediate quark with momentum p emitting a boson which is later 
absorbed by the final state quark (of momentum k). Both light front time ordered contributions to V in Eq. |2(] are 
depicted in Figure]^. Therefore, at leading order in the retarded interaction, we have shown the first step in the iterative 
scheme produces the correct potential computed from time ordered perturbation theory. 



2.2 Second order interaction 

We now illustrate the procedure at second order. At two iterations, we have the following equation for the approximate 
light front wave function 

4>W =2i?+x(l-x)z 2 | ^ ^ ^G(k,R)V(k,p)G(p 7 R)V(p, q )G(q 7 Rh(z 7 ^ rcl \M 2 ), (22) 

where z = q + / ' R + . At this point, there are a few different ways to proceed. We could first perform the q~ integral 
using our previous result Eq. (|l^) (with the necessary renaming of variables) and then the k~ integral. Carrying out 
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Sgn(Im(polc)) 


y > x 


y > x 


y < x 


y < x 






y > z 


y < z 


y> z 


y < z 


v <o 


Positive 




a b v q 




a b v k v q 




Negative 




Vk 






< y < 1 


Positive 


b 


b V q 


b v k 


b V k V q 




Negative 


a Vk v q 


a Vk 


a v q 


a 


V > i 


Positive 




<t> 


Vk 


4> 




Negative 


a b v k v q 




a b v q 





Table 3: All possible combinations of signs of the poles' imaginary parts for the second order bound state equation using 
< z < 1. The p~ poles are labeled as they appear below Eq. (p2[). A </> denotes an empty intersection of the restricted 
regions of x, y and z. 

the integration in this particular order shows one that contributions to are from the region where < z < 1. 
The resulting integrand is quite complicated and further integration and interpretation in terms of time ordered graphs 
require algebraic gymnastics. Furthermore, rather delicate cancellations are required to show that the regions y > 1 
and y < make a vanishing contribution. We shall not take this route, however, since the integrals all converge, each 
route must yield the same answer. Hence we shall use the fact that only < z < 1 contributes to perform instead the 
p~ integral first. There will be fewer combinations to enumerate and the p~ integral is essentially an example considered 
by Ligterink and Bakker who prove the equivalence of covariant perturbation theory and light-front time-ordered 
perturbation theory for scattering states. The difference is that we do not assume that contributions to the bound state 
wave function come from plus momentum fractions between zero and one, rather we derive this fact. 

There are four poles relevant to the p~ integral: two from the propagator (p~ a and p~ b ) and the remaining from the 
two iterations of the interaction (p _ l; (k _ ) and p~„(q~), see Eq. (|l2|) and make the obvious replacement of variables). 
For ease we shall abbreviate the latter as P~ Vk and p - ,^ , respectively. The signs of the imaginary parts of the poles are: 
Sgn(Im(p" Q )) = - Sgny, Sgn(Im(p" fe )) = - Sgn(y-1), Sgn(Im(p" 1 ,J) = - Sgn(y-x) and Sgn(Im(p~ 1)(f )) = - Sgn(y-z). 
Given these and < z < 1, Table || lists all possible combinations of signs depending on x,y and z. We anticipate that 
the contributions from y outside of (0, 1) will vanish. Let us then tackle these first. 

For y < 0, only one region potentially contributes to the integral: y > x and y < z. Since < z < 1, the last 
inequality is automatically satisfied. The first two inequalities, however, constrain x < 0. Closing the contour in the 
lower-half plane, we pick up only the residue at p~,, fc = p~ v (kr) which spreads the k~ dependence to more denominators. 
There are now five poles in k~: k~ a , k~f, which come from the propagator G(fc, R), two poles which come from evaluating 
the propagator G(p, R) at p~^ fc , and a final pole from evaluating V (p, q) at p - ^ . The exact form of these poles is of no 
concern, only their imaginary parts. The poles of the propagator G(k, R) have imaginary parts — — and —-^j which are 
both positive since x must be negative. The remaining poles' imaginary parts are shifted by evaluating at p _ ^ fc since it 
has an imaginary part —^37- Thus the poles which come from G(p, R) have imaginary parts — ^ 4- and — + 
which are again both positive since y < and y > x. Finally the interaction V(p, q) evaluated at p~ = p~,„ fc has a k~ 
pole with imaginary part — ^ — I — 3- and since y < z and y > x, this is positive. When integrating over k~, all poles 

y z y x 

have positive imaginary parts and thus the integral vanishes in this region. 

On the other hand, when y > 1 only one region can potentially make a contribution: y < x and y > z. Quite 
similarly to the above, since < z < 1 we must have x > 1. Closing the contour in the upper- half plane we pick up 
the same residue, namely at p ^ fc • Then in proceeding to evaluate the k~ integral, we have the same poles as before. 
This time however all the relevant signs are reversed and consequently all k~ poles lie in the lower-half plane. The 
contribution from this region is also zero and hence the only possible contributions to the second order wave function 
are truly from < y < 1. 

Referring back to Table |3], there are four distinct regions that will make a contribution to the wave function for 
< y < 1. Two are considerably simpler than the others since we can close the contour to enclose only a single pole. 
We shall proceed with these. When y > x and y > z, the pole in the upper-half plane gives us +2iri Res(p~ b ). Since 
this pole did not stem from either interaction, there is no spreading of k~ or q~ dependence to other denominators. As 
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B 





D 




E 





Figure 3: Light front time ordered diagrams contributing to the second order wave function. Notice that diagrams A-D 
can be summed into the square of the first order interaction. 



a result, the two minus momentum integrals factorize into separate integrals both involving a propagator and a single 
interaction. Ignoring the overall factors, the relevant terms are 



1 



1 



1 



(k--k- a )(k--k- 6 )(k- 



J p 6-p a (q -p b)(q -q J(q -q b) 



(23) 



Notice the factor of Dyj(y, p ± |M 2 ) in the middle corresponding to free particle propagation between interactions. The 
k~ integral is now identical to that in Eq. ( |l5| ) provided we include the factor 9(y — x). Thus < x < 1 and we pick up a 
contribution from the pole k~ a . The q~ integral is essentially the same picking up only the residue at q~ a . Taking these 
residues of Eq. ( p3| ) gives us the time ordered diagram A of Figure ||[ Parallel to the case just considered, the region 
y < x and y < z exhibits the same features. We pick up the residue of a single pole, p~ a , which does not contaminate 
the other terms with k~ or q~ dependence. The minus momentum integrals again factorize with a Weinberg propagator 
sitting between the two interactions. The light front time ordered graph B in Figure [| corresponds to this contribution. 

The last two regions of integration are more complicated for two reasons. Firstly we pick up the residues of two 
poles, one of which spreads minus momentum dependence to other energy denominators. Secondly, each residue alone 
does not correspond to a time ordered graph, rather there are two graphs entangled in their sum. Consider the case 
where y > x and y < z for which we close the contour in the upper-half plane to pick up (considering only the pieces of 
the integrand altered by evaluating at the p~ poles) 



Res(p b ) +Res(p ) oc 



1 



1 



(P 6-P „J(P 6-P a)(P 6- 
1 



1 



(P 6-P bJ(P 6-P a)(P v~P v a ) (P 6-P „J(P v -P a)(P v~P ej 



(24) 



After the algebraic manipulation that brings us to the second line, the first term has been cast in a familiar form fea- 
turing the free propagator between interactions. The k~ and q~ integrals then factorize into forms considered previously. 
Evaluating these integrals yields an expression corresponding to the time ordered diagram C in the Figure. The second 
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term, however, has one energy denominator (p~^ ~P~ v k) which has both kr and q~ dependence brought about by 
two iterations of the interaction. 

Of the two remaining integrals, consider first the one over q~. The propagator G(q, R) sitting in Eq. ( p^ ) brings along 
poles q~ a and q~ b . The first has an imaginary part — — which is negative for < z < 1, while the second has a positive 
imaginary part —zrzj- The other two poles are found in Eq. (p4[). One is of a familiar form since p - ,^ — p~ a = q~ — q~ Va 
where Sgn(Im(q _ t , a )) = — y (z~ y ) ■ This imaginary part has a negative sign since 0<z<l,0<y<l and y < z. The last 
q~ pole which stems from the denominator p - ,, — p _ ^ fc has an imaginary part + ~y~£ which is also negative. Hence 
we pick up +2ni Res(q~ b ) by closing the contour in the upper-half plane. Collecting all the terms we have organizes us 
to perform the last integration over kr 

% - x)9[y(l - y)]0[z(l - z)]9{z - y)V(qr b , z, g ± ;y, p x ) B w (z, |M 2 ) 

(k--k- a )(k--k- h )(p-Jk-)-q-„J(k--k- t ,J ■ 1 > 

Now since y > x, we must have x < 1. When x < all four k~ poles have positive imaginary parts and when < x < 1 
only k~ a has a negative imaginary part. Hence finishing up the integral requires — 2iti Res(k~ a ) and we have additionally 
shown < x < 1. The result takes the form (with 0[x(l — x)]6[y(l — y)]6[z(l — z)] understood) 

D w (x, k x \M 2 )6(y - x)V(k- a , x, k x ; y, p^V^x, k ± ;y, p ± ; z, <^)6{z - y)V{q b , z, q^; y, p x ) D w (z, q x |M 2 ), (26) 

where we have 

Having spelled out this last energy denominator which has two force carrying bosons simultaneously propagating, 
we refer the reader to Figure || to observe that diagram F corresponds to Eq. (p6|). The remaining case in Table | 
is evaluated similarly to the case just worked. There are two residues picked up and hence an algebraic manipulation 
required to untangle diagram C and diagram E. With all possible contributions to the second order wave function, we 
must now summarize the results thereby deriving the improved integral equation for the bound state. 

Each term we have evaluated has an overall factor of 2R+ p-^ - k _ ^ = Dw(i, k^rei |M 2 ) which sits outside of the p 
and q integrals. There is also an overall factor of 2 R+(q-\ _ q - ) = Dw(- Z ,q i rc i l-^ 2 ) which multiplies 7(z,q ± rcl |M 2 ) — 
these can be combined to form ip(°\z, q ± rel ). Consider the terms leading to diagrams A-D of Figure || Writing these 
out explicitly (omitting an overall Dyf(x,'k.- L ic i \M 2 )9[x(l — x)] and factor of 9[y(l — y)]6[z(l — z)] in the integrals), we 
find| 

/ 2(2%yX- v) 2(2% d z(l- 7) ( % " ^ P±) ^ P± z) * {<Ca ' *' ^ V ' P±) 

+ 6{x - y)V(h- b , y, p ± ;x, k x ) D w (y, p X \M 2 )0(z - y)V(sf b , z, q x ; y, p x ) 
+ 9{x - y)V(k- b , x, k x ; y, p x ) D w (y, p x \M 2 )9(y - z)V(q~ a , z, q x ; y, p x ) 

+ 9(y - x)V(k~ a , x, k x ; y, p X )B w (y, p x \M 2 )9(z - y)V{q~ b , z, q ± ;y, p x )^ V (0) (z, q^, (28) 

which can be rearranged as 

/ mtwtt \ (°( x - y) V{ ^^ k± ->y> p x ) + % ~ x Wk-a, x, k^ ; y , d w ( v , p x \m 2 ) 

J 2{2nyy(l - y) \ J 

dz d q / \ 

-i y 9(z-y)V(q- b7 z,q ± ;y : p ± ) + 9(y-z)V(q- a ,z,q ± ;y,p ± )j^(z,q ± ) 

''"'' pA -V(y,p x ;x : ^)^(y,p x ). (20) 



2(2tt) 3 z(1 - z) 



2(27r) 3 y(l - y) 



3 Notice the reversal in the order of the variables p and q in the interactions. This is only because we evaluated the p integral first which 
sat in the middle of Eq. Such are the joys of notation. 
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In the last line, we used the first step in the iteration scheme to relate 'ip^ to tp^ and the definition in Eq. (^0|). At this 
point, the second order wave function depends on both the first order and zeroth order wave functions. To derive the 
bound state equation, we merely change all wave functions to second order since any difference is at least third order. 
Finally, we have derived the second-order bound state equation for the light-front wave function 

where the interaction V to second order (and hence three sets of variables) is^ 

V(x, k^; y, p ± ;z, q x ) = 2(2vr) 3 2/ (l - y)5(z - y)S(q ± - p^Vfa k x ; y, p x ) 

+ 9(x - y)V(k~ b , x, k- 1 ; y, p ± )Vfe(x, k ± ;y, p x ; z, c^^y ~ z)V(q~ a , z, q- 1 ; y, p^) 

+ 9(y - x)V(k~ a , x, k ± ;y, P x )Vf(x, k x ; y, p ± ;z, q x )0(* - y)V(sC b , z, q ± ; y, p X ), (31) 



with V(x, k 4 -; y, p 4 -) as the one-boson exchange potential ( p0[ ) and Vf as appears in Eq. (p7|). The term Ve which comes 
from diagram E is 

2R+Vv{xM;y^z^)- 1 = R --^ n - ip± 2 ~ + ^ )2 _ + z f ik m+ P { 7-yf {R - k) °- (32) 

Now that we have found a scheme to derive the light front wave function from time ordered perturbation theory 
with a covariant starting place, it is time to derive the scheme in a formal setting. This will enable us more easily to 
consider matrix elements of current operators. 



3 Formal reduction 

In this section, we largely follow the reduction scheme set up in Hto convert the covariant Bethe-Salpeter equation into 
a three-dimensional light front version using the results of section |2| to motivate the choice of auxiliary Green's function. 
We must first adjust our notation to be more in line with theirs, where needed. Above, we have removed overall four 
dimensional, momentum-conserving delta functions, e.g. our propagator G(fc, R) is the momentum space version of G(R) 
defined by (R'\G\R) = (2tt) 4 ^ 4 ) (R' - R)G(R). Additionally we have removed the momentum conserving delta function 
between initial and final states of particle one in the disconnected Green's function, i .e. (k\G(R)\p) = (2ir) 4 8^ 4 \k — 
p)G(k,R). Using (R'\V\R) = (2tt) 4 J( 4 ) (R - R')V(R), we see the kernel V(k,p) defined above is V(k,p) = (k\V(R)\p). 
Notice the dependence on the total momentum R drops out since we are specializing to the ladder approximation, cf 
Eq. (|Tl|). In terms of these fully four dimensional quantities, the Lippmann-Schwinger equation for the two-particle 
transition matrix T appears as 

T = V + VGT. (33) 

A pole in the T-matrix (at some R 2 = M 2 , say) corresponds to a two-particle bound state. Investigation of the pole's 
residue gives the Bethe-Salpeter equation ([j]) for the bound state vertex T 

T = VGT. 

The Bethe-Salpeter amplitude ^ is then GT. Above, our ^Bs(k,R) is the momentum space version of with delta 
functions removed (as is similarly true of the vertex T(k, R) which we implicitly used above in going from Eq. (|l|) to Eq. 

). Following m, it is convenient to denote quantities able to be rendered in position or momentum space with bras 
and kets. We will employ this notation only for quantities that have been stripped of their overall momentum-conserving 
delta functions, for example T(k,R) = (k\Tji), where R is used as a label for the bound state for which R 2 = M 2 . 

4 This is of course the potential to second order in the ladder approximation. The full kernel contains additional second order pieces, such 
as a crossed ladder term, which contribute to the light front kernel beginning at second order. 
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3.1 The scheme 



To reduce the Lippmann-Schwinger equation to a three-dimensional light front version, we must introduce an auxiliary 
Green's function^] G in place of G (as in |l9|). To inject G, we rewrite Eq. ( |33| ) in the form 



T- 1 = V- 1 - G 



= V~ 1 -(G-G)-G, 



(34) 



ve- 



in which V 1 — (G — G) has the form of an inverse interaction (denoted W 1 ) with respect to the auxiliary Green's 
function G. Re-inverting T, we find that 



T= W + WGT, 

provided that 

W = V + V(G - G)W. 
This gives us an alternate way to arrive at the bound state vertex function 

r = WGT. 



(35) 



(36) 



(37) 



Until we choose G to reduce the Bethe-Salpeter equation, we have only formally created a more complicated problem 
to solve. 

To choose a light front reduction, G must inherently be related to projection onto the initial surface x + = 0. For 
simplicity, we denote the integration Jdk~(kr \0(R) = O(R). With this notation, we will always work in momentum 



space for which the only sensible matrix elements of 
O(R) 



0(R) arc of the form (k+,k ± \ 0(R)\ p-,p+,p ± ). The operator 

G(R) = G{R) 



is defined similarly. For a useful reduction scheme (one that preserves unitarity), we must have 

Bearing in mind an extra delta function in G(R) relative to section ||, we have already calculated G(R) in Eq. (|) 
above, thus 



(k+,k ± | 5 (i?)|p+,p ± )^(k+,k ± | G(R) |p+,p ± > = 5(k+- P +)5(k ± -p 



± ^2vri 0[x(l - x)] 
2R+x(l - x) 



DwOr,k\ el |i? 2 ). (38) 



Furthermore, we wish to choose G in order that our reduction scheme reproduces time ordered perturbation theory. Yet 
in section || we were able to recover time ordered perturbation theory without explicit recourse to an auxiliary Green's 
function. Therefore, we wish to keep G as close as possible to G. Referring back to the intuitive reduction scheme, (once 
we insert momentum conserving delta functions between the two particles) we realize that we must have tacitly used 



G{R) =G(R) 
G(R) = G(R). 

This simplest way to achieve this requires 

G(R) = G(R) g~\R) G(R), 
for which the explicit mention of G disappears. We can then write the reduced T-matrix in the form 
tW^g-^R) G{R)T(R)G(R) g~ l (R) = g^(R) G(R)T(R)G(R) g~\R), 



(39) 



(40) 



(41) 



5 Of course, the new Green's function must maintain the unitarity of the theory and thus is not completely arbitrary. 
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where t is a three-dimensional auxiliary transition matrix. As a function of the total momentum _R, if T has a bound 
state pole, so does t. This enables us to relate the four-dimensional wave function to a three-dimensional one. To this 
end, we write 



where the reduced auxiliary kernel is 



t(R) = w(R) + w(R)g{R)t(R), 



W (R) =g- 1 (R) G(R)W(R)G(R) 



(42) 



(43) 



Taking the residue of Eq. @ at R 2 = M 2 , gives a homogeneous equation for the three-dimensional vertex function 7 



\ lR ) = w(R)g(R)\j R ), 



(44) 



from which we can define the light front wave function \tp R ) = g(R)\lR)- By iterating the Lippmann-Schwinger equation 
for T twice, it is possible to relate T to t and thereby construct T given t, which is clearly not possible from the definition 
( ffl|) . Taking the residue of this relation between T and t yields the three-dimensional to four-dimensional conversion 
between bound state vertex functions, namely 



\T R )=W(R)G(R) \ lR ). 
Lastly, we can manipulate the four-dimensional Bethe-Salpeter amplitude into the form 



(45) 



**>=!+( G(R) -G(R)) W(R) G(R) \j R ) , 



(46) 



which justifies the interpretation of \ipR) as the light front wave function since 



\y R )=[g(R)+( G(R)- G(R))W(R)G(R) \\ 7r ) = \^ R ) 



(47) 



We comment that all light front reduction schemes when summed to all orders yield the x + — projection of the 
Bethe-Salpeter equation. The result of this section is stronger since the kernel is calculated in light front time ordered 
perturbation theory, e.g. w(R) calculated to second order is identical to our result of section || and carried out in this 
formal notation in ||. Lastly the normalization of the covariant and reduced wave function is discussed in Appendix A. 

3.2 Reversal of logic 



we now reverse 



Having used the intuitive reduction scheme to motivate the form of the light front reduction found in 
the logical circle and use the formal reduction to derive the iterative scheme of section |^. 

The heart of the intuitive scheme lies in integrating out the minus momentum dependence of the covariant wave 
function. So we merely rewrite this order-by-order in the formal reduction scheme. Utilizing Eqs. (|37|) and (|3^) while 
omitting total four-momentum labels since they are all identical, we can write 



1*) 



GWGlg- 1 ] 



n=0 



G 



|7>- 



(48) 



From truncating the series in n at some j < 00, we are led to a natural sequence of better approximations to the light 
front wave function provided the interaction is weak. Working out the two lowest orders to parallel the discussion above, 
we have 



\tf> 



GVG | 7 >, 



(49) 
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Figure 4: Graphical depiction of the irreducible five-point function T^ 1 as sum of impulse terms and a gauged interaction. 
By construction, T M is gauge invariant. 



which looks like one iteration of the Bethe-Salpeter equation and 



(2)\ 



GVG 



GVGVG 



GV 



It), 



VG 



It) 



(50) 



(after having used Eqs. (37) and (|45j) to get to the second line) looks like two iterations of the Bethe-Salpeter equation. 
One might wonder how the instantaneous approximation appears in any of this. Using again Eqs. ( |37| ) and (45), we 
can show 



and notice 



G 



GVG 



t) = g | 7 ) 



(51) 



It) 



(GV) 2 G | 7 ) = 



GV G |r) 
(GV) 2 G\T) 



(GV) j G | 7 ) = (GVyG\T). 

Thus at any order in the scheme of section it is as if we have iterated j times to arrive at (GVy G\T) and replaced 
the lone G to the immediate left of T with G. The instantaneous approximation then requires GIT) to have no minus 
momentum poles besides those of the propagator which is just the statement made by equation (pl|). 

With a formal reduction (that agrees with our intuitive scheme) in hand, we are now ready to tackle computing 
matrix elements of currents between bound states and related applications. In turn, these applications will help us to 
grasp better the nature of the reduction scheme. 



4 Current in the reduced formalism 

In this section, we extend the formalism presented so far to include current matrix elements between bound states. We 
do so in a gauge invariant fashion following |^o| . Our notation, however, is more in line with the elegant method of 
gauging equations presented in [pi"! . This latter method is far more general and extends to bound systems of more than 
two particles. 

Consider first the full four-point function G^ defined by 

G (4) = G + GTG. (52) 

For later use, it is important to note that the residue of G^ at the bound state pole R 2 = M 2 is n)^ Using 
the Lippmann-Schwinger equation for T, we can show the four-point function satisfies 

G {i) =G + GVG {i) . (53) 



14 



Figure 5: Graphical depiction of the five-point function G^ M . The irreducible five-point function is the gauge invariant 



To discuss electromagnetic current matrix elements, we will need the three-point funtion function d% where the label 
i denotes particle number (the utility will be clear momentarily). Now we define an irreducible three-point function T^ 1 
in the obvious way 

< = d^d,. (54) 

Now we need to relate the one-particle electromagnetic vertex function to the T matrix. Let eP = 9 M denote the 

electromagnetic coupling to our scalar particles. Since the electromagnetic three-point function is irreducible, we 
have 

Tf = G^G^df (55) 

and by using the definition of G^ Eq. |5^, we have the desired relation 

Tf = cP + TGd». (56) 

Notice the right hand side lacks the particle label i. In the first term, the bare coupling acts on the zth paticle while in 
the second term the bare coupling does not act on the ith particle.. For this reason we have dropped the label which 
will always be clear from context. 

In considering two propagating particles' interaction with a photon, the above definitions lead us to the impulse 
approximation to the current 

r o = Ti^ 1 + d r lr 2- (57) 

Additionally the photon could couple to the particles in the midst of their interacting. Define a gauged interaction V 1 
topologically by attaching a photon to the kernel in all possible places. This leads us to the irreducible electromagnetic 
vertex defined as (see Figure ^) 

= rft + v, (58) 

which is gauge invariant by construction. 

Lastly to calculate matrix elements of the current between bound states it is useful to define a reducible five-point 
function (see Figure |5|) 

G (5) M = G (4) r P G (4)^ ( 5Q ) 

Having laid out the necessary facts about electromagnetic vertex functions and gauge invariant currents, we can now 
specialize to their matrix elements between bound states. 

4.1 Form factors at leading order 

Consider Figure ||[ Knowing that the four-point function G^ (by way of the analytic structure of T) has a bound state 
pole, we can easily deduce the current between bound states. It is clear from the figure, that electromagnetic scattering 
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p P' 

Figure 6: The leading order diagram for the electromagnetic form factor 

between free two-particle states is described by the current J% = G^ 1 G^ ^G^ 1 . From this it is easy to generalize to 
scattering between bound states (labeled by P, i = initial and P',f = final) by taking the appropriate residues of the 
five-point function, namely 

J% oc Res (G (5) P 2 = M 2 , P' 2 = Afj) . 
Equality is achieved by removing \^p>) and (^p|- By using Eq. (|59|), we can cast the current into the form 

Jfi = (*p'|r M |*p>. (60) 

The vertex function satisfies the Ward-Takahashi identity, so that one can easily use the Bethe-Salpeter equation for ^> 
to demonstrate current conservation: (P 1 — P) ■ Jfi — 0. 

We now specialize to the case where the initial and final states are the same particle. Lorentz covariance and current 
conservation dictate^ 

J» = -i(P'f> + P»)F(t), (61) 

in which appears the electromagnetic form factor F(t), with t = [P' — P) 2 . In this section, we calculate the form factor to 
leading order in perturbation theory. Typically one chooses the + component of the current as well as a frame in which 
A + = (where A is defined by A M = P' M — P M ). As we plan to focus on nonvalence contributions as well as generalized 
parton distributions, we do not choose A + = which puts us in a situation similar to Sawicki p^ ], who considered such 
expressions earlier. We do, however, choose to work in a frame which simplifies the transverse momentum dependence, 
namely one in which P- 1 = 0. 

At leading order, we can remove any explicit dependence upon the interaction from the current, thus can be 
neglected. Similarly rf = d M + TGd^ ~ cP since the Born series for T starts off at V . Since our bound particles are 
equally massive scalars there is no essential difference in attaching the photon to particle one versus particle two. Hence 
at leading order we have 

J^ = 2^ pl \^d^ P ). (62) 

Instead of calculating the covariant form presented above, we use the results of section || to express the leading order 
form factor in terms of the light front wave function, via equation (|4^). But to leading order, \^) « G Using the 
plus component of the current results in 

F(t) = — i— ( 7P ,| |G(P')d + ^ 1 G(P)| | 7 p), (63) 



6 Additionally in the light front formalism, there can be an additional spurious form factor B(t). This is due to the presence of an extra 
four- vector in the problem which we introduced by needing to specify the orientation of the light front plane. For a composite scalar, 
B(t) can be removed simply by looking at the + component of the current Jf [Q. This is a manifestation of the linearity of the angular 
condition for scalar particles. 
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where we have defined £ = — A + /P + . Inserting the effective resolution of unity, we can convert the above expression 
into 

F{t) = P+(l-C/2) / W) 1 (0 I<7P ' |P+ ' P±> (P\ G ( P ') d+d 2 X G(P)\k) (k + ,k^ | 7 p). (64) 

This expression simplifies since the disconnected Green's function still possesses an overall delta function (p\G(P)\k) = 
(2tt) 4 5^ 4 ' (p — k)G(k, P). After further insertions of unity, we have the factor 

(p\d + \k) = -i(2k + +A + )(27r) 4 <5 (4) (p — k — A). (65) 

Now define x = k + /P + . With this parameterization and performing the last integral rendered trivial by the above 
delta function, the expression for the form factor reads (see Figure ^J) 

F(t)= 1—772 / (^7*(y,k^+(l- 2/ )A i |Af 2 Mfc + A)(2x-C)d(fcKP-fc) 7 (x,k^|Af 2 ), (66) 

where for the final meson's reduced vertex, we have employed y — + — jz^- Since the vertex functions are reduced, 

there is only k~ dependence in the propagators which enables integration. The poles are situated at 

k~ — k~ —is 

K a — K on ~ 

k" 6 = p- + (k- P)" - ^ (67) 
[k- c = -A- + (fc + A)- n -^ 

We cannot assume, however that contributions are only for < x < 1 and < y < 1, since the form factor appears as 
a matrix element of reduced vertices rather than light front wave functions. We went to some labor to show that the 
wave function has support only between zero and one. The reduced vertex need not vanish for x outside of zero to one, 
rather utilizing Eq. (Q) 

1 (y,p ± \M 2 ) = (yR+,p 1 -\ lR ) 

2(27032(1 -s) "^' P±; X > k± |M2) Dw(x ' k± |M2)7(:C ' k± |M2) (68) 

demonstrates that 7 for < x < 1 uniquely determines 7 everywhere. 

Evaluation of the integral over k~ in Eq. ( |66|) is straightforward (and carried out in [[l3)). We pick up only 
contributions from < x < 1, y on the other hand is between and one. There are two distinct regions of the 

integration, < x < C and £ < x < 1 (or equivalently, < y < and < y < 1 respectively) which correspond to 
picking up different residues. For x > C, the time ordered diagram which we pick up corresponds to Figure ^ (which has 
the luck of being covariant on one occasion and now time ordered). And in this case, since both x and y are between 
zero and one, both the initial and final state vertices can be described by a light front wave function. Evaluation yields 

Contribution from this region survives the limit £ — > 0. Contribution from the other region of integration vanishes in this 
limit (modulo possible zero mode singularities). In this limit, which corresponds to the usual frame for the evaluation 
of light front form factors, Eq. (^9|) reproduces the Drell-Yan formula for the electromagnetic form factor [p3[ . 

In the other region x < C, however, a Z graph confronts us in extracting the form factor, see Figure M. As in our 
previous work, we arrive at an expression 



where to define the vertex (since y < 0) we must appeal to Eq. (|68|). If interpreted graphically, we have fixed up 
the non-wave function vertex by appealing to crossing symmetry. This crossed interaction, however, involves a time 
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p P' 

Figure 7: The Z graph which confronts us to leading order in evaluating the electromagnetic form factor 



ordering with a parton traveling back through time, and ignores the spectator of momentum k. As precarious as it 
sounds, this corresponds to hadronization off the quark line and is precisely what the formal reduction dictates — at least 
from iterating the reduced Bethe-Salpeter equation. But there is obviously a reason we have not spelled out the exact 
functional form of this contribution. Here we are ascertaining the form factor at leading order. Using Eq. ( |68| ) to cross 
the interaction, introduces a power of the weak coupling and thus this nonvalence contribution is at a higher order than 
which we are currently working. [| To be correct, the region < x < £ then contributes nothing at leading order, but 
leaves a puzzle for the first order corrections. 



4.2 Form factors at next to leading order 



Above we have derived the leading order expression for the electromagnetic form factor. We were confronted with 
a nonvalence contribution, which we temporarily dealt with by iterating the reduced Bethe-Salpeter equation (fssf). 
The process of iteration, however, launched us into a higher order of perturbation theory and thus as far as leading 
order expressions are concerned, the contribution from < x < C vanishes. Here we work at next to leading order in 
perturbation theory to unravel the valence and nonvalence contributions to the form factor. 

At leading order in the weak coupling, we must retain the first term in the Born series as well as the leading term 
in converting from the Bethe-Salpeter wave function to the reduced vertex. Again there is no contribution from the 
gauged interaction as it enters at the next order. Referring to equation (pCl), we have 



J" 



2< 7 H \G(P')(l + V(P')(G(P') G(P'))) ( 



•#0 + SJ i + SJ f + 5J 1 + °\ y2 



in which appears the leading order result 



./"doT 1 + U(-A)G(-A)#d 2 - 1 j (l + (G(P) - G(P))V(P)) G(P) | 7P ) 

(70) 

as well as the first order correction terms given by 



2<7P'I 
2<7P'I 
«? = 2< 7 p,| 



SJ? 
5J» 



G{P')d"d^ L [G(P) - G{P))V{P)G(P) 
G(P r )V(P r )[G(P') - G{P'))d^d^ l G{l 
G(P')(v(-A)G(-A)dAd^G(P) | 7P 



\1p) 



\lp) 



(71) 



SJj is the correction term due to iterating 



The notation is reminiscent of the intuitive reduction scheme of section]^, e.g 

the Bethe-Salpeter equation for the final state once and then making an instantaneous vertex approximation (modulo 
two-particle reducible contributions) , while 5Jtf is the correction from retaining one order of the four-dimensional Born 
series. Indeed we could have derived the above form using the intuitive scheme. Our task is now clear: perform the 
minus momentum integrals to ascertain what contributions the new poles introduced by iteration will make. 



7 By using Eq. ( |68[ ) for the vertex, one might wonder how the interaction is small since gw <■ 
are dealing with a non-wave function vertex. Having crossed the interaction, we are describing 
reduced into the final state wave function. 



' 1 for the bound state. The key is that we 
a hadronization potential which cannot be 
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4.2.1 Born term 



Naturally, we opt to consider the simplest first. Let us write out the momentum representation of 5J^ appearing in Eq. 

(0), 

5J+ = 2i J (2^4 (2^)4 +(1 ~ v)A x \M 2 )d(j> + A)d(P - p) 

x V( P ,k)d(k)d(k + A)(-i)(2k + +A + )d(p) 7 (w,p ± \M 2 ), (72) 

where w = p + / P + and v = jEt- Part of the k~ integral should be attacked first since it vanishes. We notice the k~ 
poles k~ a , k~ c , and k^^q - ) all have negative imaginary parts for £ < w < x and hence the contribution to the form factor 
vanishes in this region. Now we can tackle the p~ integral with ease. Since the interaction is a correction to the bare 
photon vertex, this boson exchange can never be reduced into either vertex via the bound state equation. Thus we must 
have an explicit power of the weak coupling in any contributing diagrams. Were there non-wave function vertices to 
consider from the initial or final state, crossing symmetry would dictate they are at a higher order and can be neglected 
here. Hence we have both < w < 1 (from the initial vertex) and < v < 1 (from the final vertex). Combining the 
two inequalities, we have the constraint £ < w < 1. Furthermore, our above consideration of the k~ integral tells us 
that only w > x is relevant. These last two inequalities fix the signs of the p~-poles' imaginary parts and the integral 
is +2ni Res(p~ h ). Taking this residue produces the energy denominators p~ b — p~ a , p~ b — p~ c which turn the initial and 
final vertices into wave functions respectively and p~ fc — p~ v . = k~ Vb — k~. 

Now in performing the k~ integration, we note that since £ < w < 1 and w > x, we must have x < 1. Note further, for 
x < all poles {k~ a , k~ c , k _ 1 , b } have positive imaginary parts. Hence we pick up only two contributions: — 2ni Res(k~ a ), 
for < x < C; and +2iri Res(k _ „ b ), for £ < x < 1. Writing these out, we have for £ < x < 1 

—2iP + f dxd k 1 - 9(x — () dwd p 1 - 9(w — x) x — ui . i . . i . 

' -V> («,P +(1 - v)A ± ) 



[2(2tt) 3 ] 2 J x(x-() vw(l-w) g 2 

V(k- C , x, k- 1 ; w, p X ) (2x - C) V(k~ a) x, k x ;w, p X )ip(w, p ± ) (73) 

where the energy denominators are given by 

"' > - x)V(k~ a , x, k x ; w, p^)- 1 = P~ - k~ on - {P t~^ ±)2 + f ~ (P ~ p)on, (74) 



2P+(w-x) v 7 v ' ' ' r ' u " 2P+(w-x) 

w - x)V(k- c , x, k^;w, p^)- 1 = P~ + A- - (k + A)" - ^'^f \^ - (P - p)~ , (75) 



2P+(w-x) v ' v ' ' r ' v ,on 2P+(w-x) 
and for < x < ( 

+2iP + f dxdk 1 - 0{C — x) dwdp 1 - 0(w — x) i ,, > . i , ^ , , i . i , . 9 , 

ra / n V n \ - tp (v,p +(1 — v)A ) Dw g, k +oA A ) 

x (2x - C) V(k~ a , x, k x ; p X )VK P X ), (76) 

where we have defined cr = a;/^. Each of these contributions can be interpreted diagrammatically. For C < x < 1, we 
have the time ordered graph G of Figure ||, and for < x < £, we have the Z graph labeled J in Figure ^. 

4.2.2 Initial state iteration 

At the next level of difficulty, let us write out the momentum representation of SJ^ appearing in Eq. (|7l|), 

5Jt = ^ J ^^i*{yM+{i~y)^\M 2 )d{P-k)d{k + A)d{k) 

x (-i)(2k + +A + ^V(k,p)d(p)d{R - p^iw^p 1 - \M 2 ), (77) 
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Figure 8: Diagrams which contribute to the form factor at next to leading order (for x > Q). 



where again w = p + / P + ■ For ease, we have left off the subtraction of the two-particle reducible contribution (arising 
from G defined in equation^). We shall need to make use of its form below and thus explicitly the missing term is 

-2< 7 j»| G(P')d^ 1 G(P)V(P)G(P) \ 1P )=-2( 1P ,\ G(Pi)d»d^G{P) \ lP ). 



The form of this term is identical to Eq. (62), except with the crucial sign difference. Thus in evaluating SJ^, we shall 
merely throw away pieces that can be reduced back into the initial state wave function. To perform the p~ integral, we 
note that w is restricted to < w < 1 at this order. If there were contributions outside this range, we must appeal to 
crossing symmetry for the initial vertex, which necessarily introduces higher order terms. 

In carrying out the p~ integral, there are three poles: p~ a , p~ & and p~ v (kr). Given < w < 1, the first has a negative 
imaginary part, the second positive and the third depends upon the sign of x — w. Taking the relevant residues produces 
an answer we are familiar with 



dp 



_ 1 (w,p ± \M 2 )V(p,k) 

(P"-P"a)(P"-P"fc) 



(78) 



referring back to Eq. ([L8|) for the definition of V (with the necessary renaming of variables). The resulting poles of the 
k~ integrand are: k~ a , k~b, k~ c , kr Va and k _ ^ b , and there are two cases to consider buried in V depending on the sign of 
x — w. The k~ integral is nearly the same as in section |], however, there is an additional pole k~ c . 

Consider first, the term which results when x > w. In the case x > £ we do not need to tackle a Z graph contribution 
and so we proceed along this course first. With these restrictions on the plus momentum fractions, the signs of the 
poles' imaginary parts are fixed and we may perform the integration by closing the contour in the upper-half plane, 
which produces +2ni Res(k~fc). Recognizing the multiplicative factor 



P~ 



A- 



(P-fc)- -(fc + A)" 



allows us to convert the final state vertex into a wave function. The interaction dependent term is k & 
a time ordered one-boson exchange. This contribution then appears as 



(79) 



which is 



- 2iP~ l 



dxdh 



_L 



{X \^) r (y ' k± +(1 ~ y ^( 2x - c ) D w(z> k± \ m2 ) 

dwdp 1 - 9(x — w) 
2(2tt) 3 u;(1 - w) 



2(27r) 3 yx( 



ffek^.p^^p 1 ) (80) 



and will be removed by an identical term stemming from G when combined with the other one-boson exchange uncovered 
below. 

For x < w, we take up the same restriction concerning the final-state vertex x > £. The poles of the k~ integrand 
are: k~ a ,k~t>,k~ c and k - ,^. Once again, the signs of their imaginary parts are completely specified. Evaluation of the 

integral by closing the contour in the lower-half plane produces — 2iri(^ Res(k~ a ) + Res(k~ c )^j , for which any individual 
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Figure 9: The remaining diagrams (characterized by x < £) for the electromagnetic form factor at next to leading order. 

residue does not correspond to a time ordered graph. Their sum, however, can be rearranged 
1 1 



(k a -k fc )(k a — k c )(k a -k „ b ) (k c - k )(k c — k &)(k c — k U J 

1 



(k (,-k a )(k a -k „J(k c -k b) (k a -k Vb )(k c -k b )(k c -k" 



(81) 



''6 y 



into two time ordered contributions. The first term above corresponds to the other one-boson exchange then free 
propagation to interact with the photon and finally absorption into the final state; in essence, it is equation (|8^) only 
with 6{w—x). Summing these two contributions allows us to utilize the bound state equation and the resulting expression 
is identical to the leading order form factor (for x > C) Eq. j62|). Hence these contributions are removed by G appearing 
in Eq. ©. 

The second term above (|8l]), however, is not removed by G. Assembling all the pieces evaluated by residues leads to 
-2iP + f dxdk 1 ^ 9(x — C) dwdn 1 - 6{w — x) x — w ,., , , , . , , ,, 



[2(2tt) 3 ] 2 J yx(l - x) w(l - w) g 2 

x (2x - () V(k- ai x, k 1 ; w, p X )VK p x ), (82) 

which is graphically depicted by diagram H in Figure |§|. 

Lastly in dealing with 5 J? we must consider the case x < C an d come to terms with the Z graph contribution 
encountered at leading order. Before toiling away at the integrals, we note: the final state vertex will require crossing 
and thus the interaction present from iterating the Bethe-Salpeter equation for the initial state must be reduced into 
the initial state vertex, else we are at higher order, which is beyond our concern. Let us spell this out to be certain. 

To utilize the bound state equation to remove the interaction at the initial state, we need < w < 1 and x > 0. 
Performing the p~ integral yields equation (|7^) above. Not surprisingly, there are two cases to consider depending on 
the sign of w — x. When w > x, only k~ a has a negative imaginary part. Evaluating the residue gives us half of the 
one-boson exchange (k~ a — k - ,^) -1 acting on the initial state wave function. On the other hand, when x < w both k~ 
and kr c have positive imaginary parts. Taking the residue at k~ c puts the boson exchange along the qq pair's path to 
be annihilated into the photon and hence cannot be absorbed into the initial state. The remaining parts of the diagram 
require partons moving backward through time and the crossed vertex. At next to leading order, this contribution is 
thus zero. Evaluating the residue at k~{, produces the factor (k~b — k _ t , a )~ 1 which gives the other half of the one-boson 
exchange potential acting on the initial state wave function. Combining the two exchanges allows the interaction to 
be absorbed into the initial state wave function. Thus we arrive at Figure ^. But this is precisely the leading order 
contribution to the form factor for x < ( and is subsequently removed by G. This is the key observation of this work. 
In light cone time ordered perturbation theory, appealing to crossing symmetry to deal with non-wave function vertices 
is an ineffectual fix. When one works correctly to a given order, the crossed interactions are removed.^] 

8 There is, however, the special case of an instantaneous interaction which we consider in Appendix B. 



21 



4.2.3 Final state iteration 



One final contribution to the form factor at next to leading order remains, SJ^ which arises from iterating the Bethe- 
Salpeter equation for the final state. Referring back to Eq. [n], we write this contribution in momentum space as 

SJ f =21 J ( ^)4 (2^)4 ^ q± -z A± )d(q)d(P' - q)V(g, k + A) 

x d{k + A)d(P- /c)(-i)(2k++A+)d(fc)^(x,k ± ), (83) 

where z = q + /P' + and we omit the term proportional to G since we are now conditioned to spot contributions it 
removes. 

First off, we know that contributions for z outside of zero to one will require crossing of the final state vertex. 
Furthermore, in this situation no utilization of the bound state equation can remove the interaction already present. 
Thus to this order, < z < 1. Since now we deal with the final state vertex, a bit of notational adjustment is required. 
The interaction present in Eq. ( |83| ) can be written as 

-Q 2 1 

V(q,k + A)~ 



2P'+(z-y)q--q- v (k-+A-y 
and the imaginary part of the pole is — The remaining poles of the q~ integral are q~ which comes from d(q) and 

z y 

has a form we are well familiar with, and q~ b which comes from d(P + A — q). To avoid confusion, explicitly we see 

q" b = P + A + (q - (P + A))^ - (84 ) 
In this form, the q~ integral can be evaluated just like many of the integrals above 

dcf 7W-*yfofc + A) k qX _ zAX)?(K k , +A , ; x }> (85) 

(q -q J(q -q &) 

Hidden in V are two cases depending on the sign of z — y, namely 

where the primes are to keep clear the dependence on the final state momentum P' rather than the P's heretofore. It 
is useful for what follows in section to spell out the primed poles 



Kz = + q 



on 2[q+-(k+ +A+)] r lfc z(z-y) ^ 



=P- + {*-(P + A))^ - ^ 2[q+ _ (k+ +A ( }] + ie pjtjg^. 

This leaves us to evaluate the k _ integral. For y > and y > z, the poles of the integrand {k~ a , k^b, k _ c fc(, - } all 
have negative imaginary parts except for k~b. Taking the residue, produces the interaction k~b — k!~ which is half of the 
one-boson exchange potential for the final state. This contribution then will be removed by G when paired with the 
other exchange. 

Now consider the case y > and y < z for which the poles are a,b,c and v' b . In evaluating the integral, we find 
(Res(k- a ) + Res(k C )V Neither of these residues corresponds individually to a time ordered contribution. After 
algebraic manipulation into two new terms, however, one corresponds to the other boson exchange for the final state 
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and is thus removed by G when paired with the term found above. The other term gives the contribution 
-2iP+ f dxdk-t- 9(x-() dzdq 1 - 9(z - y) y-z 



312 



f dxdk 9(x-t) dzdq 6{z-y) «Z£^ (je>q x_ zA ± ) y (k - c|y>k i +A ±. Zjq x ) 

J y x \^ % ) z) q 

x (ix - C) V(k- a , y, +A^; z, q^x, k x ), 



[2(2tt) 



which can be identified as diagram I of Figure |8[ 

Finally there is one remaining case, y < (which automatically implies z > y). Only kr a has a negative imaginary 
part and hence we pick up — 27ri Res(k~ a ). The final contribution at next to leading order is 

+2iP+ f dxdk 1 ^ 0(C-x) dzdq 1 - 6{z-y) x ± 



-ip*(z,q ± -zA^ 



[2(2^) 3 ] 2 7 a(l - a)(l - x)( z(l - z) 

x Dw(a, k^ +aA ± \A 2 ) [2x-() VQf a , y, k^ +A^; g, q^x, k x ), (89) 



which corresponds to the time ordered diagram K in Figure ^|. Figures ^ and ^ then exhaust the contributions to 
the form factor at next to leading order. Thus in time ordered perturbation theory, we see there are no contributions 
from the non-wave function vertex its presence as a reducible contribution was subtracted. We comment that current 
conservation is perturbatively improved for form factors at next to leading order over the Drell-Yan formula. Lastly if 
one is solely concerned with form factors, one can obtain simpler expressions by taking £ — > 0. 



5 Application to GPD's 

Having worked through matrix elements of the electromagnetic current in a frame where A + ^ 0, we can now make 
the connection to the generalized parton distributions (GPD's) for this model. These distributions, which in some 
sense are the natural interpolators between form factors and quark distribution functions, turn up in a variety of hard 
exclusive processes, e.g. deeply virtual Compton scattering, wide-angle Compton scattering and the electroproduction 
of mesons f24| . The scattering amplitude for these processes factorizes into a convolution of a hard part (calculable 
from perturbative QCD) and a soft part which the GPD's encode. Since light cone correlations are probed in these 
hard processes, the soft physics has a simple interpretation and expression in terms of light front wave functions (25). In 
this section, we recast our results for form factors in section |] in the language of GPD's and the light cone Fock space 
expansion. Additionally one can obtain these results directly from time ordered perturbation theory using two-body 
projection operators as explicated in Appendix C. 

The GPD for our meson model is defined by a nondiagonal matrix clement of bilocal field operators 

F(XX ^ = J L(2x-0 (VtV| 1^ d+ l(°) l^>> ( 9 °) 



where q{x) denoted the quark field operator. Comparing to the current matrix element defined above (61), the definition 
of the GPD leads us immediately to the sum rule 

~ ^F(x,(,t)dx = F(t). (91) 



2-C 

Hence one can calculate these distributions from the integrand of the form factor. In this way, the light cone correlation 
defined in Eq. (^o|) has a natural description in terms of light front time ordered perturbation theory, e.g. for x > £ the 
relevant graphs contributing to the GPD are in Figure ^ and |], and those for x < ( are in Figure ||. Also note, the £ 
dependence disappears after the x integration present in the sum rule due to Lorentz invariance. We were aware of this 
in section ^ when we opted not to choose a frame in which A + = for the form factor when we knew at the end of the 
day F(t) would be a function of t = A 2 alone. 
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5.1 Continuity 

Conversion of the contributions to the form factor above into GPD's is straightforward using Eq. (|9l| ) and (63); we 



merely remove -iP+(2x ~ () and / dx from Eqs. @, @, ©, (||), @ and f dx^f from equation @. In 



2-C 

order for the deeply virtual Compton scattering amplitude to factorize into hard and soft pieces (at leading twist), the 
GPD's F(x 1 C, t) must be continuous at x = £. Continuity at the crossover is an even more pressing matter since all 
current experimental efforts to measure such functions are limited to the crossover p6f . The leading order expressions 
are continuous. This is easy to see since the contribution for x < £ is identically zero. The valence contribution is 
a convolution of wave functions one of which is ip*{y, . . .) which is probed at the end point since y — > 0. From the 
bound state equation Eq. ([l9|), we see the two-body wave function vanishes quadratically at the end points. Taking 
into account the overall weighting factor of y -1 , the valence piece vanishes linearly at the crossover. At leading order 
then, continuity is maintained at the crossover, while the derivative is discontinuous. Working only in the valence sector, 
valence quark models will never be of any use to experiment since the value at the crossover requires one wave function 
to be at an end point. In the three-body bound state problem (e.g. the nuclcon), the wave function again necessarily 
vanishes; however, the valence GPD will vanish only if the three-body interaction vanishes at the end points (which is 
physically reasonable). 

Let us now check the next to leading order contributions to the GPD for continuity. First we shall deal with the 
term stemming from iterating the Bethe-Salpeter equation for the initial state (|8^). Since there is no Z graph generated 
from iterating the initial state, we expect this contribution to vanish; after all, there is no related piece for x < £ to 
match up with. Looking at the expression, we see again ip*(y, ...)/?/ which vanishes linearly as x — > £. Moreover there 
are the interaction terms 9{w — x)V(kT a , x, k" 1 ; w, p ), which is finite as x — > £, and 

9(w - x)V(k- Cl x, ; w, p X ) X = C j-j-J^S 2 , (92) 

which vanishes at the crossover. Thus not only does the initial state iteration term vanish at the crossover, its derivative 
does so as well. 

Now we investigate the Born terms ([Z3|) and ( |76| ) at the crossover. Approaching £ from above ([73]), we encounter two 
interactions V(kr a , . . . ) which is finite as above and V(kr c , . . .) which vanishes linearly (enough to cancel the weight 
(x — C) _1 - Thus we are left with the finite contribution to the GPD at the crossover from the Born term 



dk 1 - dwdp ± 9(w-C) 



[2(2?r) 3 ] 2 J C vw(l - w) 

x r (v, +(i - ^ 1 ) (k i +A l )2 + m / ( kl 2 y . C; P X ) ^ p x ). 03) 



On the other hand, approaching the crossover from below (f76j) we have to deal with singularities as a = x/C, — > 1. 
Writing out the propagator for the quark anti-quark pair heading off to annihilation, we see 

Dw (^ + ^|A^- (k± / A - + m , (94) 

The linear vanishing of this term is enough to cancel the singular weight (1 — cr) _1 . Taking the limit x — > £ then produces 
equation ( p3|) and thus the Born terms are continuous. 

Lastly we must see how the final state iteration terms match up at the crossover. Having spelled out the Born terms, 
the final state terms follow simply once we note V(k~ a , y, k- 1 +A ± ;z, q- 1 ) is finite as y — > and 

d(z - y)V(k~ c , y, k^ +A^; z, q^) -> ^ ^ - - 2 , (95) 
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which vanishes as x goes to £. Thus at the crossover we have the contribution 



final 



[2(2tt) 



rfk J 



dzdq 1 



C(l-C) z(l-z) 



. I 2 n 

~ / k +m 



(k-L +A-L) 2 



r V 



2P+C 



,0,k^+A^,q x UP(z,^). (96) 



Approaching £ from below, we again utilize equation (|9J|) in taking the limit of (|89|). The result is ( |96| ) and hence we 
have demonstrated the GPD's continuity to first order in the weak coupling, i.e. 



F((,(;,t) = F(c,(,t) Boin + F(C,Ct) 



final 



(97) 



no matter how we approach x = C- 



5.2 Fock space representation 

We now proceed to rewrite the GDA's in terms of Fock component overlaps. In the diagonal overlap region x > £ this 
will be a mere rewriting of our results, while there is a subtlety for the nondiagonal overlaps. First let us rewrite the 
term to zeroth order in the weak coupling (p9). Define the two-body Fock component as 



1 



^XiX 2 



^(X^k^rd), 



(98) 



noting that the relative transverse momentum can be defined as k^rei = Xz\Cy — xik%. In terms of Eq. (pq), the GPD 
appears as 



lo _ 00= - 



where the primed variables are given by 



3 = 1,2 



(99) 



XjXj 



and the integration measure is given by 

JV 



[dx]x — Y[ d,Xi 5\1 — ; 



~3 - 1-C 

k'> = k} 



N 



(1-4) A± 



i=l 
N 



N 



[dk±]jv= [ 2 (2.)V- i n dk ^(^ k ' 



(100) 



(101) 



(102) 



Notice the sum over transverse momenta in the delta function is zero since our initial meson has = 0. The sum over 
j in Eq. (^9|) produces the overall factor of two for our case of equally massive (equally charged) constituents. 

To cast the next to leading order expressions for x > £ in terms of diagonal Fock space overlaps, we need to write 
out the three-body Fock component. Looking at the diagrams in Figure ^, it is constructed from the two-body wave 
function 



ip 3 (Xi,k±) 



2(2tt) 3 fx£ 



[dyWv A 



%i - x 1 )x z 8(y 2 - x 3 )6(p2 - k^-)V r (fc 1 on , x\, k±; x 3 , kg 



+ d(U2 - x 3 )xiS(yi - xi)6{pi - kj L )V r (fc 3 on , x 3 , k^; Xi,k± 



(103) 
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where i runs from one to three and the label j, which stems from the integration measure, runs from one to two. The 
awkward looking factor of \Jx2j 'g 2 arises from our definition of V. It serves to remove a factor of g since only one 
interaction is needed to place us in the three-body subspace and when combined with the prefactor in V produces 
the factor l/y/x^ which then symmetrically appears with the overall constant in -03. We discuss how to obtain this 



three-body wave function directly from time ordered perturbation theory in Appendix C. Using -03 in Eq. (103), the 
terms in the GPD at first order in the weak coupling can then be written compactly as 

F(x,c,tf h ° = / [dx] 3 [dk% J2 5 -^^rMX ± )M^M)- (104) 



j — 1,3 X/^j^^ 



One can verify that the diagrams in Figure are generated by (104). Note well that our definition of z as a momentum 



fraction with respect to the final plus momentum P' + already includes the factor of (1 — £) 1 with respect to P + 



Additionally there is a fourth diagram generated by Eq. (104) which does not appear in the figure. This missing diagram 



is characterized by the spectator quark's one-loop self interaction. Thus in performing calculations in perturbation 
theory, we should also use the one-loop result for the renormalized self interactions f(k 2 ) appearing in the spectator's 
propagator. As far as our development is concerned, we imagined this contribution tacked onto Figure ^, while it appears 
now as part of the three-to-three Fock component overlaps. The absence of this diagram does not affect continuity at 
the crossover. The missing diagram vanishes at x — C since the final state wave function is ijj*(y, . . . ). 

Now we must come to terms with the nondiagonal overlap region, x < £. At first order in the weak coupling, the 
diagrams of Figure ^ correspond to four-to-two Fock component overlaps. We have been cavalier about time ordering, 
however. The expressions we derived in Eqs. [76] and |8^ do not actually correspond to time ordered graphs. Both terms 
contain a product of time ordered propagators: one for the two quarks leading to the final state vertex and another for 
the quark-anti-quark pair heading off to annihilation. But for an interpretation in terms of a four-body wave function, 
all four particles must propagate at the same time. This is a subtle issue as a graph containing the product of two 
independently time ordered pieces (where one leads to a bound state vertex) corresponds to a sum of infinitely many 
time ordered graphs. Writing out the terms that concern us, we can manipulate as follows 

-(2P+) 2 C(1 - ODw^k 1 +aA ± \t)D w (z,q ± -z^\M 2 ) = (105) 

k c -k a q b ~q a 



1 



q a+ k c-k a I q b -q a k c -k a 



(106) 



In this form, we have produced the correct energy denominator for the instant of light front time where four particles 
are propagating. Multiplying this denominator by the three-body wave function yields the four-body wave function (up 
to constants). This is only the part of the four-body wave function relevant for GPD's (there are additional pieces for 



two-quark, two-boson states, see Appendix C). In the resulting sum (106), the first term will produce the two-body 



wave function for the final state and we will have a genuine four-to-two overlap. We do not write this out explicitly. 



The second term in Eq. (106), however, contains again the propagator for the pair heading to annihilation. Using 
the light front Bethe-Salpeter equation for the vertex (which contains infinitely many times) we can introduce a factor of 
the time ordered interaction. The resulting product of independent time orderings can again be manipulated as in Eq. 



(106). The result produces another overall four-body denominator which contributes to the four-body Fock component 



of the initial state. Since we iterated the interaction, however, this new contribution is no longer at leading order and 



can be neglected. Thus the second term in (106) does not contribute at this order. 

Having manipulated the GPD's into nondiagonal overlaps for x < £, we must wonder if continuity at the crossover is 
still maintained. In the limit x — > £ the light front energy of the struck quark goes to infinity (since it has vanishing plus 
momentum). Consequently k~ c , which contains this on shell energy, also is infinite and dominates the four-body energy 
denominator. This is identical to the reasoning above in Eqs. |7^ and ^ where instead of the four-body denominator, we 
have D\v(c, -l-erA^li) which is dominated by k~ c at the crossover. In both cases we arrive at the expressions found 
above for the crossover (^) and (|9tj). 

Having cast our expressions for generalized parton distributions in terms of the Fock components generated to first 
order in the weak coupling, we can enlarge our understanding of the sum rule and continuity at the crossover. Both 
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must deal with the relation between higher Fock components. The way C dependence disappears from (91) mandates a 
relation between the diagonal and nondiagonal Fock component overlaps that make up the GPD. The relation between 
Fock components must follow from the field theoretic equations of motion. The continuity condition itself is a special 
case of the relation between Fock components, specifically at the end points. Above we have seen our expressions are 
continuous (and non vanishing) at the crossover and explicitly that the three- and four-body components match at 
the end point (where x — C — 0). This weak binding model for behavior at the crossover is a simple example of the 
relations which exist between Fock components at the end points [p7f . More general relations must be permitted from 
the equations of motion to guarantee Lorentz covariance (e.g. in the structure of the Mellin moments of the GPD's, of 
which the sum rule is a special case). 



6 Application to GDA's 

In section |] we worked through spacelike form factors in the reduced Bethe-Salpeter formalism, then in section |5| we 
connected these results to GPD's. Below we tackle timelike form factors to demonstrate the versatility of this approach 
and make the connection to the generalized distribution amplitudes (GDA's) for this model. Analogous to GPD's, 
GDA's encode the soft physics of two-meson production and can thus be thought of as crossed versions of the GPD's. 
The GDA's enter in convolutions for various two-meson production amplitudes [ p8| . These distribution functions as 
well as timelike form factors are a theoretical challenge for light front dynamics, since there is no direct decomposition 
in terms of meson Fock components alone. Furthermore, we shall see the leading order expressions are nonvalence 
contributions (which necessarily excludes a description in terms of most constituent quark models). 
The timelike form factor F(s) for our model meson is defined by (see Figure |lC|) 

(* p VI r M \Q)=-i(p-p'rF(s), (107) 

where s = (p + p') 2 is the center of mass energy squared. Now define = p^ + p' p and £ = p + / P + . We can work 
out the kinematics of this reaction in a frame where = 

s 

2P+ 

(l-C)s 
2P+ 

,s(l-C)C-Af 2 , (108) 
where M is the meson mass. 

Similar to GPD's, the GDA for our model has a definition in terms of a nondiagonal matrix element of bifocal field 



P- = 

P~ = 
P ±2 = 
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operators 



$(Z ' C,S) = / 2tt(2z-1) Vl 9(s_)i 9+,?(0) |0> - (1 ° 9) 



Such a definition of the GDA leads directly to the sum rule for the time like form factor 

22-1 



2C-1 



${z,(,s)dz = F(s), (110) 



and hence a means to calculate $ from the integrand of the timelike form factor. 

Taking the appropriate residues of the five-point function, we arrive at Fig. ^ for the timelike form factor. Keeping 
only the leading order piece of the electromagnetic vertex r M , we have 

$( Z) C,s)=2P+| ^^ 7 *(£,k^p^M 2 )G(M^^ 

(111) 

where we have made use of z = k + / P + . Since this term carries no overall factor of the coupling, we must have only 
wave function vertices (knowing well that any apparent non-wave function vertices will be subtracted in the complete 
next order expression). This requirement translates to < | < 1 and < < 1> an d hence we do not pick up a 
contribution at zeroth order in the coupling. 

To work at first order, we pick up three terms analogous to those in section |j. We denote these as 5 J!f, 5J£ and 5J£,. 
The first Born term for the three-point electromagnetic vertex 5Jt^ is quite simple. Any diagrams contributing to this 
term must contain an overall power of the weak coupling since the boson can never be reduced into either final state 
vertex. Thus we must restrict the plus momentum fractions of the two final state vertices to be between zero and one. 
For the same reason as the zeroth order result, momentum conservation in the diagram shows this restriction cannot 
be met. This leaves us to consider only diagrams that arise from iteration of the Bethe-Salpeter equation of either final 
state meson. 

Considering first the term 5Jg, we have the contribution to the GDA 

$ p (z,C,s) = 2iP+ J 0^1^w, q ± -wp ± )G(q,p)V(q±)G(k,p)G(k~p, V ')d-\k-p) 1 *(yM-a-y)P ± ), 

(112) 

where we have chosen to abbreviate w = q + / p + , y = fE?- We have customarily omitted the subtracted term containing 



G since we know it will only remove the non-wave function vertices present in Eq. 112 which have the same for as those 
encountered at zeroth order. Thus we can entirely omit any non-wave function vertices, knowing well they will be 
canceled. Thus we have < w < 1 and £ < z < 1, which will be of incredible aid in evaluating the integrals. 

With these restrictions, the q~ integral can be performed similar to equation (|l3| ) and the subsequent k~ integration 
resembles that in section || but here all the imaginary parts' signs are fixed. Thus we get one contribution $ p to the 



GDA from 



, . . , 2 [d^e{z-Q dwdq 1 - ... ± ± . 

[2(271-)^ J yz(l — z) w(l — w) 



x v(k~ b , J,k ± ; W ,q ± )Dw(z,k ± | S )^(y,k ± -(l- 2/ )p ± ), (113) 



where for clarity we spell out the form of the interaction 



2p+U-w) V C > 2(k+-q+) 
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L M 
Figure 11: Leading order diagrams for the timelike, pion form factor. 



As far as contributions to the timelike form factor are concerned, we can interpret Eq. ( |113| ) as the time ordered diagram 
M of Figure 0. 



have 



At first order in the weak coupling, we have one term remaining to consider SJ£,. Omitting the subtraction of G, we 



$ p ,(z,C,s) = 2iP H 



dk-dk-*- d 4 q 



(2tt) 4 (2tt) 4 ' VC 



M 2 ^G(k, P )d~ 1 (p-k) 



x G(P - k, p')V(P - k, q)G(q, p')7* (w, q x +w p ± |M 2 ). (115) 



Since we are removing the crossed interactions by hand, we restrict both < w < 1 and < z < f. The signs of 
all poles' imaginary parts are now fixed and the integration proceeds just as many above. We thus arrive at the final 
contribution to the GDA 



3y(z,C,s) 



[2(2*0 



:il2 



dk^ 0(C-z) dwdq A 



er(l — cr)C w(l — w) 



r(^,k x -ap x ) 



x D w (z, \s)V(k- a , 1 - z, - k x ; w, q^iw, q x +w p x ), (116) 



where we have defined a — z/Q and the interaction is 

g 2 e((-w) 



Vik-^l-z,-^^,^)- 1 =P- -k- on -[(P-k)-q]- n -q 



(117) 



In this form we recognize this contribution as diagram L of Figure |Tl]. 

Having found the leading non- vanishing contribution to the GDA namely $ = $ p + <ly , we observe that the higher 
Fock components derived in section || (as well as in Appendix C) do not fit naturally into (112) or (116). One needs 
a Fock space expansion for the photon wave function in order to have an expression for the GDA in terms of various 
Fock component overlaps. With the expressions derived for the GDA we can use Eq. ( |110| ) to obtain the timelike form 
factor. 



7 Conclusion 

We have motivated the light front reduction of the Bethe-Salpeter equation which yields a valence wave function 
calculated from a light front, time ordered kernel. Intuitively this is performed by iterating the covariant equation any 
number of times followed by a single instantaneous approximation that allows for all of the light front energy integrals 
to be performed. Formally it is achieved by using an auxiliary Green's function G defined in Eq. (fl0|). The reduction 
scheme for the Bethe-Salpeter equation is then applied to matrix elements of the electromagnetic current. Writing 
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out the gauge invariant current completely Eq. ^8], we were able to investigate a bound state's interaction with an 
electromagnetic probe systematically in perturbation theory. We illustrated applications to form factors, generalized 
parton distributions and generalized distribution amplitudes. 

To carry out our analysis, we were forced to adopt the ladder approximation for the interaction (|ll]). With this 
model, we saw how further iteration of the covariant Bethe-Salpeter equation leads to a better approximation to the 
energy poles of the vertex function, cf Eq. for the effective two-body interaction, which has a second order term 
stemming from the four-dimensional nature of the kernel. Alternately, we can view the energy poles of the vertex as 
giving rise to contributions from the nonvalence Fock space, i.e. the additional poles which allow for the second order 
correction to the interaction also have a higher number of quanta propagating during intermediate light cone times. To 
make this connection more transparent, we saw how the energy poles generated the higher Fock sp ace contributions to 
the form factor in section |i] (written out explicitly in|^). Additionally in Appendix A, we saw in Eq. 131 how corrections 



to the normalization stemming from poles of the vertex could be interpreted as a nonvalence probability (149). 

In obtaining matrix elements of the electromagnetic current, we were confronted with Z graphs containing non-wave 
function vertices, see Eq. (||), for example. Within time ordered perturbation theory, where there are higher Fock 
components contained in the Bethe-Salpeter wave function, these non-wave function vertices are removed when one 
works completely to any given order in perturbation theory. We contrasted this with the case of an instantaneous 
interaction (or most often encountered, contact term approximation for the Bethe-Salpeter vertex) in Appendix B. In 
this special case, one can successfully appeal to crossing in order to interpret non-wave function vertices; however, in 
this approach, there are no true higher Fock components. At first mention, it would seem that one could use these types 
of models to get some glimpse of nonvalence physics. This might be the case for timelike processes, but for GPD's the 
distributions will always vanish at the crossover which is useless for comparing with experiment. Lastly one could use 
an instantaneous reduction scheme to reduce the Bethe-Salpeter equation and still appeal to crossing to have non-zero, 
non-wave function vertices. In this situation, one is sacrificing accuracy in determining the light front wave function, 
since the potential is not calculated in time ordered perturbation theory. Moreover, the higher Fock pieces then become 
entangled in this instantaneous reduction, whereas one believes that corrections to the form factor from higher Fock 
components are systematically smaller for higher particle content. From the above analysis, we expect cancelations of 
the instantaneous terms in favor of higher Fock states. 

Having worked within time ordered perturbation theory, however, our expansion of the form factor (and generalized 
parton distributions) at leading and next to leading order was economical. In this form, we were able to demonstrate that 
our GPD is continuous at the crossover (while the first derivative is discontinuous). We also unraveled the generation 
of higher Fock states from the two-body sector using our explicit expressions for GPD's. The diagonal overlaps could 
be written down directly, whereas to uncover the nondiagonal overlaps an infinite sum of time orderings had to be 
dealt with using efficacious algebra. The expressions for higher Fock contributions agreed with those calculated from 
old-fashioned perturbation theory in Appendix C. Moreover, we saw in section [| how the higher Fock states lead to the 
non- vanishing at the crossover and are hence essential for any phenomcnological modeling of these distributions. Lastly 
we investigated timelike form factors and the related generalized distribution amplitudes in this model. Expressions for 
these were systematically evaluated in perturbation theory which is in contrast to the non-existent Fock space expansion 
for these types of processes. 

With the formalism explored here, one could use phenomenological Lagrangian based models to explore both general- 
ized parton distributions and generalized distribution amplitudes within the light front framework. Of particular interest 
would be the breaking of Lorentz invariance (in particular, rotational invariance, see p9|] ) characteristic of Hamiltonian 
theories. This breaking would manifest itself as a violation of the sum rules (and polynomiality conditions). Such an 
investigation is interesting not only for testing phenomenological models, but also for anticipating problems for any 
approximate non-perturbative solutions for the light cone Fock states. Nonetheless more model studies are warranted 
before truly realistic calculations can be pursued. 

Acknowledgment 

This work was funded by the U. S. Department of Energy, grant: DE-FG03 - 97ER41014. 



30 



Appendix A: Normalization 



In section |^ the relation between the four dimensional Bethe-Salpeter wave function and the reduced light front wave 
function is presented. We did not, however, address the normalization of cither. Given the relation between reduced 
and covariant quantities, we can use the Bethe-Salpeter wave function's normalization condition to derive the reduced 
version. First let us review the covariant wave function's normalization [ p0| . 
The four-point function defined in Eq. (p3|) has the behavior 



G^(R) 



R?-M 2 +ie 



finite, 



(118) 



near the bound state pole. In what follows, we remove the total momentum labels on operators since they are all 
identical, R. Since the four-point function satisfies the equation 



G (4) = G Wfc-i_y) G W 



the Bethe-Salpeter amplitude must satisfy 



1 = lim —i- 

R 2 ^M 2 



R 2 -M 2 



(119) 



(120) 



which is necessarily finite since \^ r) satisfies the bound state equation: I^r) — GV^r). Application of l'Hopital yields 
the covariant normalization condition 



2iR» = (V R \ 



d 
dRu 



(121) 



with evaluation at R 2 = M 2 understood. 

The normalization condition for the reduced wave function is then deduced by using the conversion (|4^) in Eq. (121). 
Hence taking the plus component 



2iR+ = (j R \ G 1 + W{G-G) 







dR- 



V 



1 + (G-G)WJG 



\7R) 



(122) 



The complicated normalization condition is indicative of the higher Fock space components contained in the Bethe- 
Salpeter wave function which must now be generated from W when we deal exclusively with the reduced wave function. 
To see this explicitly, we again work in our perturbative model. 

Let us start with the leading order contribution to the reduced wave function's normalization. 



2R+ 



(ir\ \g( 



d 



OR 



■G G 



|7*> = 1- 



(123) 



To perform the integration, we note 



dR 



-G-^(k,R) = -2ii?+d- 1 (fc)(l -x), 



(124) 



where we have customarily chosen x = k + / R + . Evaluation of the integral in equation (123) is now standard given the 
two poles k~ a and k~f,. Closing the contour in the lower-half plane for x € (0, 1), we pick up the (simple) residue at k~ a 
avoiding the double pole at k~{,. Sorting out the overall factors gives 



LO 



dxdh 



2(2tt) 3 x(1 - x) 
a simple overlap of the two-body wave function. 



r(x^ ± )^(x,h ± ) = h 



(125) 
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To analyze the normalization to first order in the weak coupling, we note that our covariant interaction is a function 
of the momentum difference alone, i.e. dV/dR' 1 — 0. As in section |], we substitute W ~ V to first order. Expanding 
equation (122) to first order we find 

N LO + SN+ .. . = 1, (126) 



where N LO is the result (125) and the first order correction is 

5N = 2~RT (7ii| \ G (^ G ~ 1 )( G ~^) VG + GV ( G ~°')(m= G ~ 1 ) G \ l7fl) ' (127) 

where by Eq. (^p|), the presence of G merely subtracts the leading order result N ho . Considering just the first term in 
the above equation (and omitting the subtraction G), we have 

d * ^ (l-x)-/*(x,k ± \M 2 )d(k)d(R-k) 2 V(k,p)d{p)d(R-p)-/(y,p ± \M 2 ). (128) 



(2ir) 4 (2ir) 4 

The p~ integral is identical to that considered in deriving the bound state equation to leading order in section |^. Thus 
using (|l8| ) we have 

\ */ .li.rfij/aj/n 1^2 f d V d P 



(l-xh*(x,^\M')d(k)d(R-kY / » f ,V(k-,x,^;y,p^)(-i)TP(y,p^. (129) 

(2?r) 4 J 2(27r) J y(l - y) 

The poles confronting us are again the same as those for the kr integral of the first order bound state equation (the 
only difference of course is that k~b is now a double pole). When y > x, we pick up the residue at k~ a similar to before, 
which not surprisingly will be subtracted by the G term. 

On the other hand, for x > y, we pick up the residue at the double pole k~j, which requires us to differentiate. Using 
the product rule, we pick up the other half of the leading order result which is then subtracted by G as well as new term 

' lnl k " ' W ~ P(x, k^e(x - y )J-J-V(k- b , Xl k^;y, p^(y, p^). (130) 



2(2tt) 3 x(1 - x) 2(2tt) 3 2 /(1-?/)^ v ' 1 v y, 2R+dk 



The relevant piece of the second term in Eq. (127), is evaluated identically up to {k <-> p}. Thus we produce instead a 
factor of 

d ~ 

a P b 

To make the end result symmetrical, we merely decode the notation: V(p~ b , y, p^; x, k x ) = V(k~ a ,x,k ± ;y,p ± ). Now 
combining the two terms and their relevant 9 functions, we can rewrite the result using the explicit form of the one-boson 
exchange potential, namely 

5N - / 2^-*) v) r(x - k±) (- Tspn^v.^))^). (mi 

Thus although the covariant derivative's action on the potential vanishes, we can manipulate the correction to the 
normalization into the form of a derivative's action on the light front, time ordered potential. With this form, we can 
compare to the familiar nonvalence probability to be discussed in Appendix C (in the frame where R + = Rr = M/y/2 
and M/V2 is the eigenvalue of the light front Hamiltonian denoted p~ below). 



Appendix B: Instantaneous kernel 

We have seen in section ||, that crossing the interaction to interpret the non-wave function vertices which occur in 
evaluating form factors leads to contributions which are subtracted at the next order in perturbation theory. In our 
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model, we have seen that the correct diagrams in the reduction scheme are light front time ordered and there are no 
extra entangled contributions to deal with from hadronization off quark lines. This still leaves a provocative question. 
What about the case of QCDi + i considered in |l4j where such interaction crossing symmetry is used? 

The effective two-body interaction in the 't Hooft model is instantaneous in light cone time, and hence is independent 
of the minus momenta. Here, we specialize to the case of an instantaneous light front potential (or equivalently the 
leading piece of an instantaneous light front reduction scheme^]) to clear up the apparent disparity in dealing with the 
non-wave function vertices. 

From our comments in the text (section |^), we see that since our effective interaction is instantaneous, the fully 
covariant Bcthc-Salpeter vertex function T is light-front energy independent. In light of this observation, we can say a 
number of things: e.g. the Bethe-Salpeter wave function 'J only has "valence" poles which stem from the propagator; 
thus, in particular, no higher Fock space components can be generated by carrying out minus momentum integrals in 
expressions for form factors. This is part of the clue for dealing with nonvalence contributions directly in this specialized 
situation. We shall elaborate this for the case of form factors. Our notation is that set out in considering leading order 
expressions for form factors in section || although we mainly write down only schematic equations in this Appendix. 
We omit total momentum labels where there is no ambiguity, i.e. when they are all the same. 

With no minus momentum dependence in the vertex function, we can write the Bethe-Salpeter equation Q im- 
mediately in a form similar to Eq. (|6^) by performing the k~ integral which picks up contributions identical to Eq. 
(i 

I7>=«ffl7>. (132) 

where v is to represent the instantaneous interaction, which by nature is already three dimensional. We are able to make 
an identical statement (to one made earlier) about nonvalence contributions to the vertex, i.e. since g cx 6[x(l — x)], 7 
for < x < 1 uniquely determines 7 everywhere. 



We can define a wave function in terms of the reduced vertex in the usual way, \ip) — g\i). The reduction in (132) 
naturally shows this wave function is indeed the x + = projection of the Bethe-Salpeter equation since 

M = 1 1*) = |c|r) 

= G\ |7> = «7|7>- (133) 

In the second line, we used the fact that the reduced vertex has no minus momentum dependence and is thus already 
three dimensional. In this form it is easy to see that iteration of the Bethe-Salpeter equation (like that carried out in 
section ||) leads to nothing new; after all, there are no poles of the vertex function to better approximate. 

Now consider the electromagnetic form factor for this special case. It is defined in terms of a matrix element between 
bound states of the irreducible five-point function 

(134) 

As our concern is only with the non-wave function vertex, we can ignore the gauged instantaneous interaction's contribu- 
tion to r M . Similarly, we need only consider the bare term in the electromagnetic three-point function to be confronted 
with such a nonvalence contribution. The relevant contribution under consideration from the above equation is then 

(7/| \G(P f )d+d^G(P,)\ | 7i >, (135) 

having used |^) = G\T) = G\ I7). This form is identical to the leading order expression considered in section [| for the 
form factor. Having already evaluated this expression, we know there are two co ntrib utions, that of Figure ||for x > £ 
and the nonvalence piece in Figure We are then in the position to iterate Eq. (132) to extend the non-wave function 



vertex's definition by crossing. Above, however, we learned to treat such a procedure with care since in the process an 
additional factor of the presumably weak coupling is added. 

9 Additionally it is true in models where the vertex is described by a contact interaction. One can probe the nonvalence sector although 
there are no true higher Fock space components. 
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To convincingly derive the form of the non-wave function vertex in this instantaneous model, we first iterate the 
Bethe-Salpeter equation for the final state (any iteration for the initial state will only get reduced back into the wave 
function). Using 

(*/| = (T f \G(P f ) = {T f \G{Pf)V{Pf)G(Pj), (136) 

we can then write the current matrix element as 

(^/lff" 1 (^/)|G(F/)U(P / )G(P / )d+d 2 - 1 G( J RO| |7<>- (137) 

Now we use the fact that our interaction is instantaneous to trivially do the minus momentum loop integral added from 
the iteration of the Bethe-Salpeter equation. The above equation converts to 

(1> f \ v(P f ) \G(P f )d+d^G(Pi)\ | 7< ). (138) 

The final minus momentum integration produces two contributions. Firstly for £ < x < 1, we are able to use the bound 
state equation for the final sta te w ave function since < y < 1. Inserting (ipf\v — (7/I, we arrive back at the original 



expression in this region Eq. ( pL 3 5|) . On the other hand, when < x < £, we have y < and appealing to the bound 
state equation is of no avail. However, what remains in equation ( |l38| ) is precisely what we would have obtained by 
extending the definition of the vertex via ( |132| ). 

Therefore, in the case of an instantaneous interaction, the non-wave function vertices can be interpreted by crossing 
symmetry as is done in JT^| . In such a situation, there are no higher Fock space contributions contained in the covariant 
Bethe-Salpeter wave function. One would also have non- vanishing, non-wave function vertices to leading order in an 
instantaneous reduction scheme of the Bethe-Salpeter equation (i.e. using another form of G than considered in section 
^ to do the instantaneous reduction). Higher order terms, however, would consist of entangled instantaneous terms 
and higher Fock components. From our analysis above, we know the instantaneous pieces would cancel in favor of 
higher Fock components — making their pressence at leading order superficial. Additionally the instantaneous reduction 
of the kernel would not have the desirable feature unique to light front time ordered perturbation theory — that is, at 
higher orders one systematically adds more particles propagating at a given instant (in light cone time) which has been 
demonstrated to lead to quicker convergence ||, [3l| even over (instant) time ordered perturbation theory. Finally we 
remark that although certainly not ideal, an instantaneous interaction approximation could give some new information 
about the nonvalcncc structure of timelike form factors and timclikc decays in the light front formalism. 



Appendix C: Two-body projection operators 

The results of this paper can similarly be achieved directly from "old fashioned" time ordered perturbation theory in 
a form which utilizes projecting onto the two-body subspace of the full Fock space. As the light front Hamiltonian 
generates light cone time translation, not surprisingly time ordered perturbation theory stems from the Hamiltonian. 
For a nice, complete discussion of this formalism for our scalar model, see [[52]. In this Appendix, we show how to derive 
higher Fock space components in this formalism, thereby demonstrating the generation of higher components from the 
lowest sector we found indirectly for GPD's and form factors in section |[ 

We write the Hamiltonian as a sum of a free piece and an interacting piece which carries an explicit power of the 
weak coupling g. In an obvious notation this is 

P~ = p-+gPf. (139) 

The free term P~ is diagonal in the Fock state basis, while the interaction generally mixes components of different 
particle number (in the scalar model we consider above, the interaction is completely off-diagonal since there are no 
instantaneous terms). Let us suppose that in the full Fock basis, we have an eigenstate of the light-front Hamiltonian, 
i.e. 

(p-+gP I -)\Tp) = p-\i>), (140) 
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where the eigenvalue is labeled by p~. Since the coupling is presumed small, the mixing of Fock components with a 
large number of particles will be small. Thus one imagines our bound state will be dominated by the two-body Fock 
component. 

To make this observation formal, we define projections operators on the Fock space V and Q in the usual sense (i.e. 
V 2 = V, Q 2 = Q, V + Q = 1, etc.). The operator V projects out only the two-particle subspace of the full Fock space 
and hence Q projects out the compliment. Let us define the action of these operators on our eigenstate 

V\1>) = \ih) (141) 

ei^> = l^e>- (1 42 ) 

Given V + Q = 1, we have immediately that l^) = \tp2) + | V'q) - We now wish to derive an equation for the two-body 
state in terms of effective two-body operators. First we note 

VP-\il>)=p-V\il>)=p-\ih), (143) 

and then write V P~l = VP^(V+Q) = Pf, v V +P^qQ, where we have defined the following notation for any 
operators A and B, P^ B = AP~B. Using this operator relation in Eq. ( |143| ) yields 

Pvq\^q)= (p--Pvv)\^)- (144) 
Now we use the same procedure on the state QP~ \ip) to find 

P Q v\^) = (p-Pqq) IV>0>- (145) 



Combining Eqs. (144) and (|145| ), we arrive at the following equation for the two-body Fock component 

P£\ih)=p-\ih), (146) 

where the effective two-body Hamiltonian is 

P- S = P vv + V cS (147) 
= p- p +P~ Q — — — Pq V . (148) 



The effective two-body interaction V c s defined in equation (147) is dependent upon the energy p~ since we have suppressed 
the degrees of freedom of the Q subspace. This yields a relation between the Q-space probability (i.e. the nonvalence 
contribution) and the effective interaction 

(■>Pq\^q) = (HV^IH ■ (149) 

In a weak binding limit, we can series expand the effective interaction in powers of the coupling and thereby re-derive 
the results of section |] for the light front reduced potential. Given that every boson emitted must be absorbed in the 
two quark sector, we can have only an even number of interactions and hence 

oo / _ \ 2n 

V ta = VgPfQ- 1 p= y £(-^=) QgPYP ■ (150) 
P r ° „=o \P r ° J 

So, for example, at leading order we have all possible ways to propagate from the two-body sector and back with only 
two interactions in between. The diagrams in Figure ^ correspond to the two possibilities distinguished by the action of 
1 p _ between interactions. At the next order, we have all possible ways to propagate from two bodies to two bodies 
with four interactions in between (see Figure |3|) , etc. 
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To generate higher Fock components from the two-body sector, we necessarily must look at the Q space state, 



I^Q) = _ p _ Pqv 1^2) • (151) 

To generate an n-body Fock component from this state, we merely act with an n-body projection operator which we 
shall denote Q n . Similar to the above, we expand in powers of the coupling to find 

P ^° n =0 \P r ° J 

For example, the leading order three-body state is obtained by attaching a boson to a quark line in the only two possible 
ways (and adding the light front energy denominator at the end). With these three body states, we can consider all 
possible three-to-three overlaps that would contribute to the form factor. These are depicted in Figure || (with the 
exception of a quark self- interaction) . The four-body sector is richer since there are two boson, two quark states as well 
as four quark states. The two-to-four overlaps required for GPD's must have four quarks. At leading order, we generate 
the diagrams encountered above in Figure ^. Not surprisingly directly applying time ordered perturbation theory from 
a light front Hamiltonian agrees with our derivation above from covariant perturbation theory in the Bethe-Salpeter 
formalism. 
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